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FOREWORD 


This  document  is  the  final  technical  report  on  a study  of 
structural  motion  effects  on  pilot  opinion  and  performance  during 
tracking  in  a large  flexible  vehicle  in  turbulence.  The  study  was 
conducted  at  the  Los  Angeles  Aircraft  Division,  Rockwell  Inter- 
national Corporation,  International  Airport,  Los  Angeles,  Calif. 
90009,  under  USAF  Contract  F33615-71-C-1418.  The  contract  work 
was  sponsored  by  the  Air  Force  Flight  Dynamics  Laboratory  under 
Project  8219,  Task  No.  821904,  and  was  under  the  direction  of 
Mr.  Frank  L.  George  (AFFDL/FGC)  of  the  Air  Force  Flight  Dynamics 
Laboratory,  Wright-Patterson  AFB,  Ohio.  The  inclusive  dates  of 
research  were  6 December  1971  to  19  August  1973. 

The  research  was  performed  under  the  supervision  of 
L.  U.  Nardi  at  the  Los  Angeles  Aircraft  Division.  Project  per- 
sonnel included  C.  A.  Crother,  Principal  Investigator; 

B.  Gabelman;  and  D.  Langton. 

This  report  was  submitted  on  19  August  1973. 

This  technical  report  has  been  reviewed  and  is  approved. 
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ABSTRACT 


Pilot  performance  parameters,  such  as  pilot  ratings,  tracking  errors, 
and  pilot  response  characteristics,  are  determined  for  two  longitudinal 
tracking  tasks,  using  a large,  flexible  bomber  in  a turbulent  environment. 

The  two  tasks,  terrain  following  and  air  refueling,  were  accomplished  on  a 
limited,  six-degrees-of-freedom  motion  simulator.  The  effect  of  structural 
motion  on  pilot  performance  and  opinion  and  the  correlation  with  pilot  ratings 
are  of  concern  for  potential  application  in  pilot  rating  prediction  methods. 
The  study  results  indicate  that  the  motion  effects  of  turbulence  and  struc- 
tural flexibility  do  not  affect  pilot  performance  or  pilot  opinion  in  the 
two  tasks  evaluated. 
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SECTION  I 


INTRODUCTION 


The  analytical  prediction  of  pilot  ratings  was  successfully  demonstrated 
for  the  hover  task  through  the  work  of  R.  0.  Anderson  of  the  Air  Force  Flight 
Dynamics  Laboratory'.  In  his  study  (reference  1),  Anderson  utilized  the  large 
mass  of  experimental  hover  data  to  develop  a relationship  between  pilot  rating 
and  measurable  pilot  task  parameters.  The  potential  significance  of  this 
methodology  is  in  its  application  to  pilot/vehicle/ analyses  and  evalua- 
tion in  terms  of  handling  qualities  specifications.  The  extension  of  this 
concept  to  conventional  vehicle  tasks  in  the  longitudinal  axis  was  undertaken 
recently  by  Anderson  (reference  2)  and  extended  to  the  lateral/directional 
axis  by  E.  D.  Onstott  of  Northrop  (reference  3).  A degree  of  success  was 
achieved,  but  specific  criteria  have  not  been  developed  to  date,  and  part  of 
the  difficulty  seems  to  lie  in  the  poor  correlation  between  pilot  ratings  and 
pilot  task  parameters  when  the  vehicle  is  subject  to  higher  levels  of  turbu- 
lence or  exhibits  poor  damping  characteristics.  Additionally,  there  is  some 
concern  that  structural  mode  motions  adversely  affect  flying  qualities  opin- 
ions. In  order  to  develop  this  methodology  for  application  to  the  large  flex- 
ible vehicle  where  structural  mode  motions  are  significant,  a study  was  made 
to  determine  the  effects  on  pilot  opinion  and  performance  of  structural  mode 
excitation  and  to  incorporate  these  effects  into  the  methodology  of  predicting 
flying  qualities  in  a turbulent  environment.  The  successful  development  of 
pilot  rating  predictions  based  on  this  data  base  would  enhance  the  potential 
for  ultimately  specifying  or  evaluating  flying  qualities  criteria  through 
these  methods. 

This  report  details  the  insults  of  that  study  in  structural  mode  effects 
on  pilot  opinion  and  performance  during  two  tracking  tasks.  Section  III 
describes  briefly  the  background  effort  in  the  prediction  of  pilot  ratings 
based  on  Anderson's  method.  Section  TV  covers  the  program  objectives,  a 
description  of  the  study  vehicle,  and  the  tracking  tasks  and  data  analysis 
methods  employed.  The  simulation  results  of  the  terrain- following  and  air- 
refueling tasks  are  presented  in  sections  V and  VI,  respectively.  Application 
of  the  simulation  results  to  the  paper  pilot  method  is  covered  in  section  VII, 
while  the  study  conclusions  are  given  in  section  VIII. 
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SECTION  II 


SENARY 


Considerable  effort  has  been  and  is  being  expended  to  improve  and  extend 
current  flying  qualities  specifications.  One  of  the  areas  of  technical  activ- 
ity has  been  in  the  use  of  prediction  methods  of  pilot  ratings  based  on 
closed- loop  analyses  and  pilot  modeling.  Limitations  in  the  method  have 
appeared  in  recent  work  for  conditions  involving  high  turbulence  with  or  with- 
out poor  aircraft  short  period  damping.  This  study  examines  the  application 
of  predictive  pilot  ratings  to  a large,  flexible  bomber  in  a turbulent  environ- 
ment where  structural  motion  is  significant. 

Two  longitudinal  mission  tasks  were  evaluated  using  a six-degrees-of- 
freedom  simulation  of  the  aircraft  equations  of  motion.  The  piloting  tasks 
were  conducted  on  a limited  six-degrees -of- freedom  motion  simulator  to  gener- 
ate a data  base  for  deriving  a suitable  pilot  rating  algorithm.  Pilot  comments 
and  performance  in  the  terrain -following  and  air-refueling  tasks,  however, 
revealed  little  correlation  between  rating  and  turbulence  or  structural  flex- 
ibility effects.  Pilot  model  characteristics  were  well  defined  for  the 
terrain-following  task  and,  in  the  resulting  lead-lag  transfer  function,  the 
lag  time  constant  correlated  well  with  turbulence  level.  In  the  air-refueling 
task,  considerable  variation  occurred  in  the  pilot  frequency  response  charac- 
teristics. A frequency  dependent,  large  time  deiay  was  found  necessary  for 
good  phase  matching  between  pilot  model  and  bode  p3ot.  Although  seme  evidence 
exists  of  a correlation  between  pilot  rating  and  pilot  pitch  stick  spectral 
density  function,  the  general  conclusion  of  the  simulation  effort  is  that  the 
pilots  rated  the  vehicle- task  problem  and  were  not  influenced  by  turbulence 
or  structural  mode  excitation.  This  is  attributable  to  the  fact  that  the  pri- 
mary tracking  tasks  were  in  terrain  following  and  ail  refueling  and  that  tur- 
bulence and  structural  effects  were  additive  to  these  tasks  and  were  not  of 
sufficient  magnitude  to  dominate  the  pilot’s  display  tracking  capabilities. 

The  current  concept  of  predicting  pilot  ratings  which  isolates  flying 
qualities  from  ride  qualities  appears  to  be  valid  based  on  these  results.  The 
parametric  variations  conducted  in  the  study  proved  to  be  of  little  influence 
in  pilot  rating  determination  and,  consequently,  the  data  base  generated  could 
not  be  used  to  develop  flying  qualities  specifications.  Unfortunately,  the 
simulation  and  simulator  availability  did  not  permit  a timely  and  beneficial 
reorientation  of  the  study  program. 
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SECTION  III 


BACKGROUND  OF  ANALYTICAL  METHOD 


The  VTOL  vehicle  exhibits  basic  instability  in  the  hover  mode,  necessi- 
tating the  use  of  stability  augmentation  systems.  Since  the  degree  of  arti- 
ficial stabilization  drastically  affects  the  vehicle's  response  characteris- 
tics, a simple  statement  of  flying  qualities  requirements  in  the  conventional 
terms  was  not  possible.  In  an  effort  to  circumvent  this  problem, 

R.  0.  Anderson  developed  analytical  methods  of  evaluation  and  prediction  of 
pilot  ratings  based  upon  closed  loop  analyses  of  the  hover  task.  Using  experi- 
mental hover  data  available,  he  showed  a correlation  between  pilot  ratings  and 
pilot  task  parameters  and  then  developed  an  analytical  method  for  predicting 
such  ratings  (reference  1).  The  developed  relationship  between  pilot  rating 
and  pilot  task  parameters  for  a hover  situation  was  found  to  be: 


R = + + 1.0 


where 


cr  - <r 


R-i  = 


m 


m 


, <r  = 0.80  = required  performance 


cr  - cr  + 10  cr 
X q 


0 < R S 2.50 


and 


= standard  deviation  of  x displacement  in  feet 


cr^  = standard  deviation  of  pitch  rate  in  rad  per  second 
R.  = 2.5  T , R S 3.25 

1 Le 
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R-  - 1.0  T.  , R < 1.20 
3 u 3 

x 

T = Pilot  lead  time  constant  in  pitch  in  seconds 

L0 

T^  = Pilot  lead  time  constant  in  displacement  in  seconds 
x 

In  an  attempt  to  expand  the  concept  to  provide  predictions,  Anderson  said 
that,  if  the  pilot  naturally  works  to  minimize  his  rating,  the  aforementioned 
pilot  rating  function  could  be  used  as  an  algorithm  for  minimizing  pilot 
ratings  based  on  the  trade-off  between  workload  and  performance.  He  applied 
this  prediction  method  to  additional  hover  data  and  obtained  good  results.  In 
an  effort  to  automate  the  application  of  the  prediction  method  and  thus  elimi- 
nate the  laborious  computational  chore,  J.  D.  Dillow  developed  the  "paper 
pilot"  which  is  a digital  computer  program  for  predicting  pilot  ratings  for 
the  hover  task  (reference  4) . 

G.  L.  Teper  of  STI  undertook  to  further  assess  the  Anderson  prediction 
method  for  the  hover  task  (reference  5) . Although  the  pilot  rating  criteria 
in  the  hover  task  used  in  Teper' s work  was  considerably  different  from 
Anderson's,  Teper  showed  that  the  Anderson  pilot  model  can  be  used  successfully 
if  the  time  delay  is  assumed  to  be  variable.  Additionally,  in  order  to  improve 
the  performance  prediction,  remnant  injection  was  found  necessary. 

Extension  of  the  paper  pilot  to  conventional  vehicles  and  a pitch  track- 
ing task  was  examined  by  Anderson  in  reference  2.  In  general,  the  results 
proved  promising,  but  the  data  base  was  admittedly  limited.  The  form  of  the 
pilot  model  and  the  rating  algorithm  used  are  as  follows: 


Pilot  Model 


where 


PR  = 


R1  + R2  + 1.0,  PERF  + R2  + 1.0  < 10 
10  , PERF  + R2  + 1.0  >10 
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0.1 

0. 974-0- 


-2.5  Tl>  Tl  < 0 


R1  = 


2.5  T , 0 £ Tl  < 1.3 
3.25,  Tl>1.3 


0.1 

0.974-  or 


, cr<  0.9739 


PERF  = 


1011  (0.974 -a)2  - 3xl07  (0.974-b)  + 3xl03,  a 2 0.9739 


One  interesting  aspect  of  the  study  results  5s  that,  for  a high  dynamic 
pressure  (high  aircraft  longitudinal  frequency)  flight  condition,  an  open- 
loop  prediction  algorithm  based  on  damping  provided  better  correlation  with 
the  experimental  data,  while  for  a low  dynamic  pressure  (low  frequency)  case 
(power  approach),  the  closed-loop  paper  pilot  did  better.  Additionally,  the 
open- loop  prediction  was  better  even  for  the  low  dynamic  pressure  condition 
whenever  poor  damping  existed. 

In  reference  3,  E.  D.  Onstott  and  E.  P.  Salmon  of  Northrop  Corp,  attemp- 
ted to  further  extend  analytical  prediction  methods  to  the  lateral-directional 
control  of  the  conventional  air  vehicle  in  a turbulent  environment.  Although 
their  work  showed  relatively  good  prediction  of  performance  measures  in  roll 
and  heading  control,  no  real  success  was  achieved  in  predicting  pilot  ratings. 
However,  the  data  base  generated  and  the  performance  prediction  methods  used 
provide  necessary  ingredients  for  further  work  in  this  area.  In  an  extension 
of  this  effort,  Onstott  et  al  further  applied  their  performance  prediction 
methods  to  two  other  fighter- type  vehicles  with  the  objective  of  generating 
flying  qualities  specifications  (AFFDL  Contract  F33615-71-C-1067) . A pitch 
tracking  task  was  also  added  to  the  study.  Results  similar  to  those  presented 
in  reference  3 wore  obtained,  and  flying  qualities  specifications  were  pro- 
posed in  terms  of  performance  measures.  One  interesting  conclusion  of  the 
study  was  that  pilot  ratings  alone  are  insufficient  for  flying  qualities  in 
turbulence  evaluations,  a position  comparable  to  that  expressed  in  reference  5 
which  indicated  correlation  with  pilot  rating  but  not  with  performance  is 
insufficient  for  satisfactory  specifications. 

Drawing  on  this  background  in  the  use  of  the  Anderson  method,  an  inves- 
tigation was  made  to  determine  the  paper  pilot's  applicability  to  flying 
qualities  specification  development  based  on  a large,  flexible  bomber  in 
turbulence.  The  results  of  that  study  are  presented  hereafter. 


SECTION  IV 


PROGRAM  OBJECTIVES  AND  SCOPE 


PROGRAM  PLAN 


The  study  objectives,  as  indicated  earlier,  are  to  define  flying  qualities 
in  turbulence  criteria  through  the  application  of  Anderson's  paper  pilot  con- 
cept. The  vehicle  being  evaluated  is  a large,  swingwing  bomber.  Two  flight 
conditions  representative  of  mission  tasks  requiring  precise  vehicle  attitude 
control  were  selected  for  the  application  of  the  analytical  methods. 

The  sequence  of  steps  which  were  to  be  followed  in  the  study  are  as 
follows : 

1.  Incorporate  the  bomber  characteristics  and  flight  controls  into 
the  pitch  paper  pilot. 

2.  Modify  the  pilot  model  and  rating  algorithm  as  necessary  to 
reflect  the  piloting  tasks  involved. 

3.  Identify  flying  qualities  criteria  based  on  these  paper  pilot 
applications  to  the  vehicle  in  turbulence. 

4.  Develop  flight  controls  for  the  vehicle  based  on  the  flying 
qualities  criteria  previously  identified. 

5.  Test  the  method  by  conducting  a simulation  program  to  verify 
the  flight  controls  design  and  the  pilot  model  and  performance. 

Unfortunately,  the  uniqueness  of  the  terrain- following  task  and  display 
necessitated  the  running  of  the  simulation  first  in  order  to  identify  the 
pilot  model  to  be  used  in  the  paper  pilot  program.  In  addition,  the  other 
task  (air  refueling)  was  sufficiently  different  as  to  require  still  another 
pilot  model.  Thus,  the  program  was  essentially  reversed,  with  the  simulation 
being  used  to  generate  data  for  the  analytical  program  and  the  flying  quali- 
ties specifications  being  based  or.  the  paper  pilot  and  simulation  results. 


DESCRIPTION  OF  VEHICLE,  TASKS,  AND  PILOT  MODEL  IDENTIFICATION 
VEHICLE  DESCRIPTION 

The  vehicle  selected  for  the  study  is  a large,  swingwing  bomber.  Its 
primary  mission  is  a low-altitude,  high-speed  condition,  and  it  uses  a struc- 
tural model  control  system  to  minimize  its  structural  flexibility  effects. 


9 


Preceding  page  blank 


The  command  stability  augmentation  systems  (SAS)  are  conventional,  and  the 
longitudinal  control  system  is  shown  in  block  diagram  form  in  figure  1.  The 
vehicle  characteristics  are  listed  in  table  I for  the  two  flight  conditions 
of  interest,  and  the  vehicle  longitudinal  equations  of  motion  used  in  the 
analytical  evaluations  are  given  in  figure  2.  These  equations,  including  the 
structural  mode  effects,  are  described  in  reference  6. 


TASK  DESCRIPTIONS 

The  two  vehicle  mission  tasks  selected,  which  require  precise  attitude 
control,  are  terrain  following  (TF)  and  air  refueling  (AR).  Since  turbulence 
effects  are  of  significance  at  these  conditions,  they  are  ideally  suited  for 
the  study.  A brief  description  of  each  follows. 

The  TF  task  occurs  at  the  low- altitude,  high-speed  flight  condition.  The 
TF  algorithm  developed  for  the  vehicle  is  the  ADLAT  system.  For  manual  con- 
trol, the  display  of  g error  to  the  pilot  is  augmented  by  a lagged  pitch  stick 
displacement  signal.  The  g error  is  generated  by  multiplying  the  flight  path- 
error  by  12.7  g per  rad.  Extensive  simulation  studies  conducted  in  support 
of  a bomber  development  program  have  shown  that  such  display  lead  information 
is  necessary  for  satisfactory  manual  terrain  following.  For  this  study,  the 
TF  task  was  restricted  to  the  pitch  axis  (vertical  plane)  with  the  longitudi- 
nal SAS  operative. 

The  simulation  tests  used  the  B-l  motion  simulator  (figure  3),  and  the 
vehicle  six-degrees-of-freedom  rigid  body  equations  of  motion  (described  in 
reference  6)  were  simulated  in  the  hybrid  computer  facility.  The  simulator 
pilots  flew  TF  tests  using  a conventional  ADI  where  the  horizontal  flight 
director  needle  indicated  the  g error  (scaled  for  1 inch  per  g).  Although 
six  degrees  of  freedom  were  simulated,  the  vehicle  was  well  behaved  lateral - 
directionally  such  that  the  pilot  task  was  restricted  to  the  vertical  plane 
only  (the  complete  lateral-directional  stability  augmentation  was  always 
operative  and  pilot  comments  indicated  no  coupling  occurred) . The  turbulence 
inputs  were  also  restricted  to  Wg  effects  only.  The  TF  task  was  flown  over  a 
specific  terrain  profile,  and  the  runs  were  from  3 to  7 minutes  in  duration. 
The  runs  always  began  at  the  same  point  along  the  ground  track  (the  statis- 
tical characteristics  of  the  task  are  described  later) . The  turbulence  was 
generated  through  filtered  white  noise  to  approximate  the  Dryden  form  of  the 
power  spectra  of  the  turbulence  and  then  recorded  on  FM  tape,  thus  providing 
a constant  turbulence  time  history,  desirable  because  of  the  relatively  short 
simulation  runs.  The  turbulence  characteristics  used  are  shown  in  figure  4. 

In  order  to  assess  the  effects  of  turbulence  or  pilot- induced  structural 
mode  excitation  on  the  pilot,  the  vehicle  simulation,  which  was  for  a rigid 
vehicle,  had  to  be  slightly  modified.  In  keeping  the  simulation  changes  as 


\ 
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TABLE  I.  VEHICLE  CHARACTERISTICS 


TERRAIN- FOLLOWING  (TF)  DATA 


Stability  Derivatives 

X = -0.02S05 
u 

Z = -0.9498 
w 

M*  = -9,006 
0 

X = -0.05644 
wg 

X = -0.0S644 
w 

Zg  = -165.79 

Mq  = -0.6517 

M = -0.00158 
Wg 

v° 

M = -0.0003 
u 

M = -0.00023 
w 

Z = -0.1055 
u 

M = -0.00158 
w 

Modal  Data 

Mode  1 

Mode  2 

Mode  3 

co  (rad/sec) 

13.59 

14.12 

21.20 

(ft/sec2) 

'i 

8.749 

-91.135 

-4.7821 

2 

Mr,  (rad/sec  ) 
i 

-0.2028 

-0.09777 

0.21949 

Er,  (1/ft  sec) 

lw„ 

2.529 

-0.9948 

-0.5867 

® 2 
Er,  (1/rad- sec  ) 

-2227.06 

-217.05 

614.67 

x8 

Ey,  (1/rad- sec) 

-136.87 

23.152 

50.266 

iq 

Er,  (1/sec2) 

7.826 

-16.673 

66.807 

\ 

E^  (1/sec2) 

13.877 

-305.86 

6.9015 

\ 

Er,  (1/sec2) 

7.755 

34.891 

-8.239 

< 

-0.0887 

0.0531 

-0.1097 

4>’g  (1/ft) 

-0.0036 

-0.0079 

-0.0030 

7). 

'l 

,»P 

(i/ft) 

0.0289 

0.0287 

0.03585 

7). 

S 

0.02 

0.02 

0.02 

TABLE  I.  VEHICLE  CHARACTERISTICS  (CONT) 


Pitch  Control  System  Gains 

K * K(h  ) = -0.4176  sec  K • K(MJ  • K(h  ) = 0.0649  rad/in. 
q p e 0'  v p' 

Knz  • K(h  ) = -0.04716  rad/g  Km  = 0.01099  1/in.,  = -0.00627  in./fps 


AIR-REFUELING  (AR)  DATA 


Stability  Derivatives 


X = -0.0049 

z = 

= -0.6121  M 

«.=  -37.188 

X = 0.02237 

u 

w 

0 

Wg 

X = 0.02237 

Z* 

= -50.44  M 

= -0.4259 

M = -0.003768 

w 

0 

q 

Wg 

X*  = 0 

M 

* 0.000042  M-  = -0.000253 

8 

u 

w 

Z = -0.0998 

M 

» -0.003768 

u 

w 

Modal  Data 

Mode  1 

Mode  3 

c o (rad/sec) 

9.0138 

12.2792 

18.7314 

2 

Z^  (ft/sec  ) 
i 

-15.805 

1.2688 

4.724 

? 

M-n  (rad/sec“) 
'i 

-0.002227 

-0.001889 

-0.01544 

Eyj  (1/ft-sec) 

g 2 
E^  (1/rad- sec  ) 

0 

0 

0 

-16.381 

-548.69 

-0.7316 

h 

Eu  (1/rad-sec) 

3.750 

-25.376 

-0.05473 

lq 

E„  (1/sec2) 

-29.732 

1.2584 

1.1284 

\ 

E,,  (1/sec2) 

\ 

-10.279 

L 

3.3782 

-27.690 

\ 
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TABLE  I.  VEHICLE  CHARACTERISTICS  (CONCL) 
Modal  Data  (Continued) 


(1/sec2) 

-0.03545 

-0.00160 

-0.04398 

A 

\ 

0.3626 

-0.0620 

-0.1500 

(1/ft) 

\ 

-0.01239 

-0.003238 

-0.01105 

i 

t — \ 

rt 

0.0198 

0.0288 

0.0363 

h. 

0.02 

0.02 

0.02 

i 

Pitch  Control  System  Gains 

= -0.6118  sec  Kg  • K(MQ)  • K(h  ) = 0.1005  rad/ in. 

) = -0.06118  rad/g  Km  = 0.0.1536  1/in.,  = -0.00087  in./fps 


- 32.20  + X y u + X^yW  + X58  + Xwgw9 

Uo  + Z(jU  + 

3 

ZyyW  + Z58  + E Zrj.rj-, 

i at  1 t 

(^Mq  + U0M\^ 

q + (mu+  MwZu)u  + ( 

r 

Mw  + Mwz 

+ (ms  + M 

wzs)s  + z rj.{ 

\ n 1 

Mw)  ^ + IV 

(*>U\  + + E^q  + 

' £E>?i„  v 
j « ’ 1 

Ucq  - w 

- E^'i 

i = i 

* q - E< 

i = 1 

G . 

Figure  2* 

Vehicle  Longitudinal  Equations 

simple  as  possible,  the  following  was  adopted.  A single,  2 Hz  structural  mode 
was  simulated,  capable  of  being  excited  by  either  gust  or  pitch  stick  displace- 
ments or  both.  The  eftect  of  the  mode  on  pilot  acceleration  was  fed  into  the 
motion  system.  The  sensitivity  of  the  mode  to  the  gust  and  pitch  stick  was 
varied  in  order  to  assess  their  effects  on  the  pilot  rating  of  and  performance 
in  the  task.  Since  the  objective  was  to  assess  pilot  motion  effects  as  a 
function  of  magnitude  of  structural  oscillations,  it  was  assumed  that  one  mode, 
with  excitation  sensitivity  as  a variable,  would  suffice.  A block  diagram  of 
this  portion  of  the  simulation  setup  is  shown  in  figure  4. 

Two  subjects  were  used  in  the  TF  simulation: 

• Subject  A has  16  years  in  military  and  civil  aviation  as  a pilot 
and  engineering/experimental  test  pilot  on  a variety  of  military 
and  civil  aircraft  and  helicopters.  He  is  presently  assigned  as 
project  pilot  for  the  B-l  program.  His  past  flight  experience 
includes  the  F-8,  SR- 71,  F-101,  F-100,  F-84F,  T-38,  B-57,  and  the 

Rockwell  Sabreliner.  He  is  a graduate  of  the  USAF  Experimental 
Test  Pilot  School  at  Edwards  AFB,  Calif. 

• Subject  B is  a member  of  the  technical  staff  at  the  B-l  Division 
of  North  American  Aerospace  Group  and,  as  an  engineer  in  the 
Flight  Control  Analysis  Group,  has  had  over  10  years  of  experi- 
ence in  evaluation  of  flying  qualities  of  fighters,  bombers,  and 
transports  simulated  in  the  Flight  Simulation  Laboratory.  He 
has  a commercial  pilot's  rating,  both  single  and  multiengine, 
with  over  1,000  flight  hours  experience. 

The  subjects  were  asked  to  evaluate  the  TF  task  as  presented  to  them  and  not 
the  motion  environment  in  which  they  found  themselves.  Hence,  their  Cooper- 
Harper  ratings  should  vary  due  to  motion  effects  only  if  their  performance 
is  impaired. 

Those  simulation  parameters  important  to  the  paper  pilot  program  and  to 
the  task  performance  measurements  were  recorded  on  FM  tape  for  later  data 
analysis.  Unfortunately,  one  sensor  (the  accelerometer  measuring  actual  pilot 
g on  the  motion  simulator)  operated  intermittently  at  best,  and  no  ronfident 
record  of  that  parameter  v/as  obtained. 
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The  AR  task  was  conducted  at  0.7  mach  at  a medium  altitude.  Unlike  the 
low-altitude,  high-speed  condition,  the  wings  are  swept  forward  and,  hence, 
the  vehicle  is  more  sensitive  to  turbulence.  To  provide  a repeatable  tracking 
task  for  the  evaluating  pilots,  a simulation  run  was  recorded  during  which 
altitude  changes  were  made  approximating  the  workload  involved  during  air 
refueling.  The  resulting  flight  path  angle  (7)  time  history  was  then  used 
as  a tracking  command  for  the  simulator  data  runs  (the  7 command  characteris- 
tics are  given  later).  Initial  attempts  at  running  the  simulation  using  the 
ADI  horizontal  flight  director  needle  as  the  flight  path  error  indicator 
proved  unsuccessful  because  of  the  conflict  in  cues  displayed  by  the  pitch 
attitude  ball  directly  behind  the  needle.  To  separate  the  task  display  from 
the  ADI , a cathode-ray  tube  (CRT)  display  was  used  with  only  the  flight  path 
error  and  cross  bar  reference  mark  being  shown  (scaled  for  1 inch  per 
10  degrees  of  flight  path  error).  Since  the  simulation  included  six  degrees 
of  freedom,  the  pilot  had  to  hold  speed  and  wings  level,  but  pilot  comments 
indicated  this  was  not  a factor  in  their  performance.  A picture  of  the  CRT 
display  is  given  in  figure  5 (the  runway  is  a TV  picture  being  provided  as  a 
background  for  demonstration  purposes) . Both  the  tracking  task  and  the  tur- 
bulence time  histories  were  recorded  and  repeated  for  the  different  runs  to 
hold  these  parameters  constant.  The  turbulence  inputs  and  structural  mode 
representations  were  similar  to  those  used  in  the  TF  task.  The  one  difference 
is  that  the  AR  data  include  turbulence  effects  alone  on  the  rigid  aircraft, 
in  addition  to  the  structural  mode  excitation  effects. 

Two  subjects  were  used  in  the  AR  simulation: 

Subject  B is  the  same  one  who  participated  in  the  TF  simulation. 

• Subject  C is  a former  engineering  test  pilot  with  27  years  of 
flight  experience.  With  a total  of  10,000  flight  hours  (5,000 
test  hours) , he  has  flown  a wide  spectrum  of  large  aircraft 
including  the  B-29,  B 50,  3-47,  B-52,  KC-135,  and  Boeing  707. 

He  has  considerable  B-52  AR  flight  experience.  As  was  done  in 
the  TF  simulation,  both  subjects  were  asked  to  rate  the  vehicle 
during  the  AR  task. 

In  both  tasks,  the  basic  problem  presented  the  pilot  is  one  of  7 tracking. 
The  pitch  stick  to  7 transfer  functions  (SAS  on)  for  the  two  flight  conditions 
are  given  in  table  II  for  the  rigid  body  aerodynamics.  The  inclusion  of  the 
three  structural  modes  would  add  three  poorly  damped  second  orders  to  the 
transfer  functions  with  the  associated  frequencies  close  to  those  shown  in 
table  I.  As  will  be  seen  later  in  the  discussion  on  the  paper  pilot  results, 
the  addition  of  the  structural  modes  does  not  significantly  affect  pilot 
motion  relative  to  the  task. 
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PILOT  MODEL  IDENTIFICATION  METHOD 


The  simulation  data  generated  had  to  be  analyzed  to  provide  task  perform- 
ance measures  and  pilot  models  suitable  for  the  respective  tasks.  The  primary 
tool  used  to  accomplish  this  is  a data  processing  and  frequency  analysis  pro- 
gram already  developed  and  in  use  at  the  North  American  Aerospace  Group.  This 
program  calculates  the  autocorrelation  and  crosscorrelation  functions  for  one 
pair  of  time  measurement  records  and  estimates  spectral  density  functions  from 
these  calculations.  Coherence  estimates  (similar  to  the  linear  correlation 
function  described  in  reference  7)  are  also  calculated  and  can  be  used  to 
identify  pilot  remnant  effects.  The  program  also  provides  the  gain  and  phase 
relationships  between  the  two  input  records  based  on  the  appropriate  spectral 
density,  cospectral  density,  and  quadrature  spectral  density  functions. 
Finally,  the  program  may  be  used  to  obtain  a complex  curve  fit  based  on  an 
iterative  least  squares  approximation  to  the  spectral  data.  The  mathematical 
expression  for  this  curve  represents  the  system  (pilot  model)  linear  transfer 
function.  The  two  time  histories  may  also  be  shifted  in  time  relative  to  each 
other  to  remove  the  nonlinear  phase  effects  of  a constant  time  delay.  A more 
detailed  description  of  this  analysis  program  is  given  in  appendix  I. 

For  use  in  the  aforementioned  computer  program,  the  simulation  time  his- 
tories recorded  on  FM  tape  were  digitized  at  50  samples /second  with  a resolu- 
tion of  about  0.05  percent.  To  prevent  frequency  folding  effects,  the  FM  data 
were  filtered  through  a 0.0072-second  filter  because  of  noise  presence  (pri- 
marily 400  Hz) . 


SECTION  V 
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TERRAIN- FOLLOWING  DATA  EVALUATION 


SIMULATION  DATA  ANALYSIS 

The  terrain- following  (TF)  simulation  data  were  analyzed  through  the  data 
processing  and  frequency  analysis  program  to  identify  those  test  characteris- 
tics important  to  the  paper  pilot  program.  A typical  set  of  data  produced  by 
the  analysis  program  is  given  here  as  an  example. 

A strip  chart  trace  of  a typical  terrain  and  vehicle  flight  path  time 
history  is  shown  in  figure  6.  Tlie  difference  between  the  command  flight  path 
and  the  actual  flight  path,  augmented  by  a signal  which  is  a function  of  stick 
displacement,  was  converted  to  an  acceleration  error  and  presented  to  the  pilot 
on  the  ADI  as  described  earlier.  Typical  time  histories  of  the  display  signal 
and  the  pilot's  pitch  stick  for  that  TF  segment  are  shown  in  figure  7.  In 
addition,  the  commanded  flight  path  angle  time  history  is  also  shown  in  fig- 
ure 8.  The  results  of  the  data  analysis  for  this  run  are  shown  in  figures  9, 
10,  11,  and  12.  Figure  9 shows  the  autocorrelation  functions  for  the  input- 
output  signals  (the  horizontal  bar  and  pitch  stick,  respectively).  Of  par- 
ticular interest  in  these  plots  are  the  autocorrelation  values  at  zero  lag 
number  since  they  are  the  variances  (p-2)  for  the  two  signals  (0.103  inch2  and 
0.017  g2).  Figure  10  depicts  the  spectral  density  functions  for  the  two  sig- 
nals, while  figure  11  shows  the  cospectra,  quadrature  spectra,  and  coherence 
functions.  (Refer  to  appendix  I for  a description  of  these  functions.) 

Finally,  figure  12  presents  the  bode  plots  for  the  input-output  signals.  This 
frequency  response,  computed  from  the  spectral  data  for  this  run,  represents 
the  linear  pilot  model.  The  two  data  points  on  the  phase  plot,  seemingly 
separate  from  the  others,  actually  should  be  plotted  at  -189.4  and  -187.1 
degrees  but,  because  the  plotting  routine  restricts  the  data  points  to  between 
±180  degrees,  such  points  are  folded  over.  The  smooth  lines  on  the  amplitude 
and  phase  plots  represent  the  results  of  an  automatic  routine  for  a curve  fit 
to  those  data  points  which  exhibit  a coherence  value  greater  than  0.5.  The 
match  is  good  except  for  the  dc  level  of  the  gain.  Th  -omputer-provided 
transfer  function  is 


3.87 


(-0.297S  * 1) 
1.38S  + 1 


while  a better  dc  gain  fit  would  be  3.4  rather  than  3.87. 


Preceding  page  blank 
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Figure  8.  Gamma  Command  (y  ) Time  History  (No.  0418) 


The  negative  lead  indicated  in  the  computed  pilot  model  is  a common 
characteristic  of  all  the  TF  runs  analyzed.  Although  at  first  glance  it  seems 
to  be  at  odds  with  most  literature  on  pilot  models,  this  model  is  not  unrea- 
sonable when  it  is  realized  that  this  is  a linear  model  being  fitted  to  data 
which  probably  exhibit  a time  delay;  i.e.,  a linear  model  including  a first- 
order  Pade'  approximation  for  the  time  delay  would  also  have  a negative  lead 
term.  As  will  be  seen  in  section  VII,  the  computed  paper  pilot  model  (first- 
order  lag  with  first-order  Pade'  approximation  for  a time  delay)  compared  well 
with  the  measured  pilot  model  of  the  form  previously  shown.  Since  the  time 
series  analysis  program  could  time-shift  the  output  relative  to  the  input, 
some  early  attempts  were  made  to  account  for  pilot  time  delays  before  the 
spectral  analyses  occurred.  However,  since  the  process  is  an  iterative  one 
(to  find  the  proper  time  delay)  and,  since  the  program  did  well  in  model  com- 
putation without  time  shifting,  further  attempts  were  ceased. 


Although  the  spectral  analysis  program  does  not  specifically  compute  the 
remnant  effects  of  the  pilot,  it  does  compute  the  coherence  function  (R)  which 
is  an  indication  of  the  linear  relationship  between  input  and  output.  As 
mentioned  earlier,  this  is  the  same  as  the  linear  correlation  function  described 
in  reference  7.  Since,  as  shown  in  appendix  I,  the  coherence  function  (R)  is 
given  as 


R2(w)  = 1 • 


Snn« 

S»(j 


where  S^(j  <j)  is  the  spectral  density  function  of  the  pilot's  total  output, 
then  j o>) , the  spectral  density  function  of  the  remnant  (noise) , can  be 
computed.  If  the  remnant  is  small,  then  it  can  be  ignored  and  the  linear 
model  justified.  If  the  remnant  is  large,  then  quite  possibly  nonlinear  or 
uncorrelated  effects  are  occurring  and  the  linear  model  is  restricted  in 
applicability.  In  general,  all  of  die  TF  runs  showed  good  correlation  as 
evidenced  by  the  coherence  function  plots  shown  in  figure  11  and  in  appen- 
dix II.  Figure  13  shows  the  PSD  of  the  computed  remnant  for  the  pilot  transfer 
function  previously  shown.  The  ratio  of  the  remnant  standard  deviation  to  that 
of  the  pilot's  total  output  (cr  snn/°sxx3  is  0.48.  It  should  be  further  noted 
that  the  time  series  analysis  program  curve  fit  only  uses  frequency  data  points 
which  exhibit  a specified  level  of  coherence.  For  the  TF  runs,  a value  of 
0.5  was  used. 


The  spectral  density  and  autocorrelation  functions  for  the  corresponding 
7C  for  this  particular  run  are  shown  in  figure  14. 
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Figure  13.  Spectral  Density  Function  of  Pilot  Remnant  (No.  0418) 


The  flight  path  command  time  history  naturally  depends  upon  how  well  the  I 

pilot  follows  the  terrain.  However,  it  can  be  expected  that  the  actual  varia-  | 

tion  in  the  statistical  description  of  this  signal  will  be  small  if  the  pilot's 
capabilities  are  not  widely  disparate  and  if  the  pilot's  performance  is  not 
significantly  different  between  runs.  Figure  15  shows  the  change  in  spectral  j 

density  function  for  one  pilot  as  the  stick  effect  on  structural  mode  excita-  a 

tion  progresses  from  nill  to  very  heavy  (curves  1 to  3).  The  corresponding  • 

autocorrelation  functions  are  also  shown  in  figure  15  with  the  resulting  | 

standard  deviations  ( <r)  of  0.065  rad,  0.077  rad,  and  0.078  rad,  respectively.  j 

Figure  16  shows  they,,  time  histories  for  these  runs.  These  variations  are  f 

typical  of  the  extremes  encountered  and,  for  application  in  the  paper  pilot  f 

program,  the  <r  of  the  Tc  used  was  obtained  by  averaging  over  a number  of  runs  | 

with  the  average  being  0.069  rad.  i 

\ 

i 

An  interesting  application  of  the  pilot  mod.  1 generated  in  this  program  | 

provided  tentative  validation  of  the  model  characteristics.  Computer  runs  \ 

were  made  with  a TF  system  using  the  pilot  model  to  close  the  g control  loop.  4 

The  run  was  made  for  an  encounter  with  an  isolated  peak.  Comparative  computer  J 

runs  were  made  using  an  automatic  TF  system.  Figure  59  in  appendix  II  shows  T 

the  two  runs  and  illustrates  how  well  the  pilot  model  performed.  * 


SIMULATION  RESULTS  ? 

The  IT  simulation  data  were  examined  to  determine  the  effects  of  turbu-  - 

lence  and  structural  mode  excitation  on  pilot  ratings  and  performance  param- 
eters. Additionally,  data  evaluation  was  made  to  identify  any  correlation 
existing  between  pilot  ratings  and  performance  parameters.  These  results  are  ! 

presented  herein.  ? 

Pilot  rating  influence  due  to  the  variables  of  gust  and  stick  excitation 
of  pilot  motion  differed  between  pilots  A and  B,  as  shown  in  figure  17.  Pilot  ! 

A maintained  constant  rating  of  2 across  the  board  of  parameter  variations.  , 

Even  for  the  case  of  maximum  effects  of  gust  and  stick  simultaneously,  his 
comment  was  "motion  doesn't  seem  to  affect  any  stick  motion  or  control."  > 

Pilot  A's  rating  of  2 for  the  task  continued  even  when  a subjective  rating  of  j 

E for  vibration  intensity  was  given  (on  a scale  of  A to  F for  increasing  inten-  i 

sity.)  Pilot  B showed  a rating  degradation  for  high  stick  effects  on  pilot  ; 

motion,  but  not  for  gust  excitations  even  though  his  subjective  rating  for  - 

vibration  intensity  was  higher  (D  versus  B)  for  the  worst  gust  case  compared  \ 

to  the  worst  stick  case.  The  effect  of  large  stick  excitation  on  performance  ] 

parameters  is  not  pronounced  as  will  be  seen  later.  The  explanation  seems  to 
tie  in  with  the  fact  that  the  stick  input  into  the  structural  mode  was  not 
rate  limited  in  the  simulation;  hence,  stick-effect  gains  large  enough  to  ' 

provide  noticeable  pilot  motion  for  small,  tracking- type  stick  displacements 
could  cause  a relatively  severe  jolt  if  the  stick  were  pulsed  or  moved 
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stepwise.  This  presumed  mental  inhibition  seems  to  have  caused  the  pilot 
rating  to  deteriorate  even  though  his  performance  did  not.  (See  figure  19.) 

Pilot  rating  is  also  plotted  versus  four  pilot  performance  parameters  as 
shown  in  figures  18  and  19.  In  the  latter  plot,  the  ratings  as  a function 
of  the  standard  deviation  ( c r)  of  horizontal  bar  (needle)  and  pitch  stick  are 
depicted.  No  obvious  relationships  exist  (note  pilot  B shows  ratings  of  2, 

3.5,  and  4 for  essentially  the  same  value  of  error  signal, ^Hg) • In  figure  18, 
the  ratings  are  plotted  versus  measured  pilot  model  lag  and  lead  time  constants. 
(It  should  be  noted  that  pilot  model  computations  were  not  made  for  all  pilot 
B test  cases,  but  only  those  to  show  the  range  of  gust  and  stick  effects.)  It 
would  be  expected  that  the  pilot  rating  would  increase  with  decreasing  rLAG 
and  increasing  tLEAD>  but  only  in  the  case  for  pilot  B is  this  seen,  even  to 
a slight  extent,  with  decreasing  rLAG 5 i.e.»  rating  increases  above  2 only 
when  rLAG  2oes  below  0.90  second. 

Performance  parameters  descriptive  of  the  pilot  model  and  task  activity 
are  shown  in  figures  20  through  24.  The  first  figure  shows  how  the  standard 
deviation  of  pitch  stick  activity  varies  with  gust  and  stick  effects  on  struc- 
tural motion.  Pilot  A shows  a slight  decrease  in  stick  °pg  with  an  increase 
in  stick  effect  on  motion.  Since,  in  this  simulation,  the  increase  in  stick 
effect  does  not  result  in  an  increase  in  vehicle  moment  control  but  only  in 
pilot  motion  through  the  structural  mode  excitation,  it  might  be  expected  that 
the  stick  activity  should  remain  relatively  constant  since  the  tracking  task 
doesn't  really  change.  (The  pilot  motion  was  never  large  enough  to  cause  con- 
sistent involuntary  stick  motion.)  Pilot  B shows  a more  constant  activity 
with  only  a slight  increase  for  both  the  stick  and  gust  effects.  However, 
pilot  A does  show  a definite  trend  in  crpg  with  gust  effects.  This  tendency 
to  decrease  stick  activity  as  high-frequency  motion  effects  increase  is  not 
unexpected.  .The  pilot  will  tend  to  decrease  the  frequency  of  his  inputs  and 
follow  the  lower  frequency  tracking  errors.  This  is  seen  in  figure  25  which 
shows  the  spectral  density  functions  for  the  three  conditions  being  discussed. 
The  curves  labeled  1 to  3 (for  increasing  gust  excitation)  shows  the  frequency 
corresponding  to  peak  power  drops  off  from  0.5  to  0.3  and  finally  below 
0.1  rad  per  second. 

In  figure  21,  the  standard  deviation  in  the  horizontal  bar  (needle)  move- 
ments is  presented  versus  the  pilot  acceleration  resulting  from  gust  and  stick 
excitation  of  the  structural  modes.  Between  the  two  pilots,  the  most  obvious 
difference  is  in  level  of  the  signal,  particularly  at  the  zero  level  for 
structural  motion.  Since  stick  displacement  feeds  into  the  display,  it  could 
be  expected  higher  stick  activity  would  result  in  higher  display  activity. 
However,  figure  20  does  not  support  that  contention  since  both  pilots  show 
about  the  same  crpg.  Examination  of  the  stick  and  horizontal  bar  PSDs  for  both 
pilots  reveals  a significant  difference,  however.  Pilot  A shows  the  same  peak 
frequency  in  both  stick  and  horizontal  bar  (figure  49)  PSD's,  while  pilot  B 


39 


I 


APILOT  A 
OPILOT  B 


STRUCTURAL  EXCITATION 

.UNSHADED  - STICK  EFFECTS. 
SHADED  - GUST  EFFECT 
TAGGED  - NO  EXCITATION 


PILOT 

RATING 


STANDARD  DEVIATION  OF  HORIZONTAL  BAR  ~ g RMS 


PILOT 

RATING 


n ~ STANDARD  DEVIATION  OF  PITCH  STICK  ~ INCH  RMS 
PS 


Figure  19.  Terrain-Following  Pilot  Rating  Versus  <7^  and  tfpg 
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(figure  54)  indicates  the  influence  of  stick  on  horizontal  bar  is  not  nearly 
as  significant.  Hence,  pilot  A's  large  cthb  can  be  attributed  to  the  large 
power  content  in  pitch  stick  at  about  0.5  rad  per  second  frequency.  This 
effect  of  frequency  is  further  borne  out  by  noting  that  the  convergence  in 
tHB  values  for  large  stick  effects  is  duplicated  by  convergence  in  PSD  simi- 
larities as  seen  in  figures  51  and  56. 

The  variation  in  pilot  model  parameters  are  shown  in  the  next  three  fig- 
ures. The  lead  time  constant  remains  relatively  steady  for  both  pilots  as 
seen  in  figure  22.  Both  pilots  converge  on  a value  of  0.30  second  with 
increasing  pilot  motion  resulting  from  increased  structural  effects.  Pilot 
lag  time  constant  shows  a definite  trend  with  increased  turbulent  motion. 

(See  figure  23.)  Both  pilots  decrease  their  lag  time  constants  from  about 
2.5  to  3.0  seconds  to  about  1 to  1.5  seconds.  Finally,  the  pilot  model  low- 
frequency  gain  (dc)  is  shown  in  figure  24.  The  gains,  normalized  to  the  no 
structural  motion  case,  decrease  with  increasing  pilot  motion  due  to  increasing 
structural  effects.  Percentagewise,  pilot  B shows  a larger  decrease  than  does 
pilot  A.  This  is  due  primarily  to  the  fact  that  pilot  B's  initial  gain  is 
higher  than  pilot  A's  (refer  to  table  IV  in  appendix  II)  because  of  pilot  B's 
greater  activity  in  the  low-frequency  range  as  mentioned  earlier. 

The  TF  simulation  revealed  the  ability  of  the  pilots  to  maintain  fairly 
constant  performance  over  the  range  of  structural  motion  effects  studied  with- 
out any  significant  trend  in  pilot  workload.  This  is  reflected  in  the  pilots’ 
unchanging  rating  of  the  TF  task  even  when  subjected  to  a high  level  of  tur- 
bulent type  of  motion.  This  is  undoubtedly  due  to  a large  extent  to  the  fact 
that  the  tracking  task  was  a low-frequency  one  and  that  the  high-frequency 
disturbance  did  not  appear  in  the  display,  but  only  affected  his  "ride." 

Hence,  in  this  respect,  the  TF  task  was  unlike  other  tracking  tasks  described 
in  earlier  referenced  works  where  the  turbulence  provided  the  tracking  problem. 

The  form  of  the  measured  pilot  model  was  consistent  throughout  the  runs, 
and  the  data  exhibited  good  coherence  throughout  the  frequency  range  sampled. 
The  negative  lead  character  ->tic  identified  in  the  pilot  model  is  not  unexpected 
if  the  assumption  is  uie  that  it  represents  part  of  the  first-order  linear 
approximation  of  the  time  delay  existing  in  the  pilot's  response. 
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SIMULATION  DATA  ANALYSIS 

Hie  air  refueling  (-AR)  simulation  data  were  analysed  in  a manner 
similar  to  tliat  done  for  the  terrain- following  (TF)  data.  Aside  from  the 
change  in  flight  condition,  the  significant  differences  in  task  were: 

1.  The  display  did  not  have  lagged  pitch  stick  displacement  summing 
Wlt1'  TbRR0R‘ 

2.  The  7 tracking  task  was  generated  in  an  open- loop  fashion. 

A typical  set  of  AR  traces  and  data  are  presented  herein  with  a discussion 
of  the  AR  simulation  data  in  later  paragraphs. 

Figure  26  shows  the  tracking  task,  y , as  a function  of  time 
[.or  - 0.88  degree  RMS  and  the  mean  = -0.^9  degree)  and  its  spectral 

density  function.  For  a typical  AR  run,  time  histories  of  the  display 
signal  fy,  RR0(J  aid  pitch  stick  are  shown  in  figure  27.  Figures  28  and 
29  show  the  autocorrelation  and  crosscorrelation  functions,  respectively. 

'Ihe  spectral  density  functions  are  presented  in  figure  30,  while  figure  31 
shows  the  cospectra,  quadrature  spectra,  aid  coherence  functions.  And 
finally,  figure  32  presents  the  bode  plots  for  these  input-output  signals. 

It  is  at  this  point  that  the  data  processing  aid  frequency  analysis  results 
differ  from  those  obtained  in  the  TF  task.  Automatic  curve  fitting  to 
these  data  was  not  achieved.  'Ihe  bode  plot  indicates  the  traisfer  func- 
tions consists  primarily  of  lead  terms.  'Ihe  analysis  progran  fa- led  to  find 
acceptable  curve  fits  for  these  lead-only  situations.  Hence,  limited 
curve  fitting  was  done  by  haul  and,  as  will  be  seen  later,  even  linear 
lead-only  transfer  functions  (as  suggested  in  this  bode  plot)  were  often 
inadequate.  Complexities  arose  primarily  because  of  the  nonlinear  phase 
behavior  (evident  in  the  referenced  figure)  of  many  of  the  plots. 


SI MULATION  RhSULTS 

As  was  evident  in  the  TF  task,  the  pilot  ratings  obtained  were  not  a 
function  of  turbulence  or  structural  mode  excitation  but,  for  the  AR  task, 
reflected  whether  the  stability  augmentation  system  (SAS)  was  off  or  on. 
Figure  33  shows  that  pilot  B rated  the  AR  task  with  SAS  on  as  a 3 regard- 
less of  turbulence  while,  for  SAS  off,  his  rating  was  a 5 for  no 
turbulence  and  a 9 for  3 fps  RMS  of  turbulence.  This  latter  rating  must 
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Figure  28.  Autocorrelation  Functions  of  Pitch  Stick  and 

Display  (No.  01A11) 
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Figure  31.  Cospectra,  Quadrature  Spectra,  and  Coherence  Functions 
of  Pitch  Stick  and  Display  (No.  01A11) 
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Figure  33.  Pilot  Ratings  Versus  Gust  Velocity  for  Air-Refueling  Task 
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be  considered  something  of  an  anomaly  because,  in  subsequent  runs  with 
SAS  off  and  with  more  extreme  turbulence,  the  task  was  rated  a S.  Pilot 
comment  and  the  data  indicate  a PIO  condition  existed  (see  figure  78  in 
appendix  III)  of  sufficient  severity  to  cause  him  to  give  a 9 rating. 

Such  a PIO  condition,  however,  did  not  reoccur.  The  time  histories  in 
the  figure  shew  the  oscillations  at  about  0.3  Hz,  and  the  pitch  stick 
spectral  density  function  in  the  same  figure  shews  the  significant  pewer 
contained  in  the  0.3  Hz  or  2 raa-per- second  frequency  region  (compare  this 
plot  with  that  shown  in  figure  79  for  SAS  off,  the  same  turbulence  level 
and  with  structural  mode  effects  but  with  no  power  peak  unrelated  to  the 
task).  Although  the  rating  seems  extreme  for  a PIO  condition,  the  pilot 
felt  he  could  not  accomplish  the  AR  task  safely.  Based  on  his  performance 
and  ratings  for  similar  conditions,  it  is  felt  this  case  is  an  isolated 
one  where,  probably,  pilot  motivation  was  temporarily  insufficient. 

Pilot  C was  likewise  consistent  in  his  ratings,  a 2 for  SAS  on  and  a 
4 for  SAS  off  regardless  of  the  level  of  turbulence  or  mode  of  excitation. 

These  data  indicate  that  the  primary  parameters  of  variation  used  in 
the  study  are  not  those  which  influence  pilot  ratings  of  the  task. 
Turbulence  level,  structural  mode -excitation,  pilot  acceleration,  and 
other  parameters  which  may  correlate  with  these  are  insufficient  to 
develop  a pilot  rating  prediction  algoritlim.  In  the  AR  task,  only  SAS 
affected  pilot  ratings.  The  gain  of  the  command  SAS  had  to  be  reduced  by 
a factor  of  4 at  the  beginning  of  the  simulation  to  provide  an  acceptable 
pilot  rating,  and  this  gain  effect  is  shewn  in  figure  34.  Increasing  the 
gain  caused  an  increase  in  stick  activity,  and  pilot  acceleration, 

all  for  zero  turbulence.  Unfortunately,  this  correlation  is  not  easily 
verified  in  the  paper  pilot  scheme  since  that  algorithm  was  developed  on 
the  assumption  that  the  optimum  pilot  input  sensitivity  had  been  selected 
and,  as  a consequence,  pilot  gain  is  not  a variable  in  the  rating 
algorithm  or  cost  function.  Examination  of  the  pitch  stick  spectral 
density  plots  for  the  three  SAS  pilot  input  gains  indicates  that,  as  the 
frequency  content  of  the  activity  decreased,  the  ratings  improved 
(figures  39,  76,  and  75  correspond  to  the  normalized  electrical  gains  of 
1,  2,  and  4 shown  in  figure  34).  The  pilot  obviously  prefers  to  work  at  a 
lower  frequency.  This  certainly  is  in  an  indication  of  why  the  TF  and  AR 
tasks  are  well  rated  by  the  pilots,  in  spite  of  their  motion  environment, 
as  long  as  the  frequency  content  of  the  task  remains  low.  The  effect  of 
SAS  on  the  frequency  content  of  pilot  activity  is  illustrated  in  figure  35 
for  runs  No.  01A11  and  04A1  which  shew  SAS  on  and  off,  respectively. 

Even  though  the  pitch  stick  <r  is  higher  for  the  SAS-on  case  (0.192  versus 
0.173  in  RMS),  there  is  a well- pronounced  frequency  peak  at  2 rad  per 
second  for  the  SAS-off  situation. 

Another  SAS  control  system  parameter  was  evaluated,  the  stick  pre- 
filter which,  however,  proved  to  have  little  effect  on  pilot  ratings  in  the 
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lag  range  evaluated.  Performance  parameters  resulting  from  the  lag 
variations  of  from  0.10  to  1.0  second  are  shown  in  figure  36.  Pilot  B shows 
relatively  little  difference  in  stick  activity  with  turbulence  or  stick  lag 
variation,  while  the  tracking  error  cr^  does  show  the  effects  of  increased 
stick  lag.  This,  however,  did  not  significantly  influence  pilot  rating. 

Additional  effects  of  SAS  on  task  performance  parameters  are  given  in 
figure  37  which  presents  the  standard  deviations  in  pitch  stick  and 
horizontal  bar  (7  ) as  a function  of  turbulence  level,  turbulence 

excitation  of  structural  mode,  and  SAS  status.  No  obvious  correlations 
exist  except  for  the  higher  o-'s  for  the  high  turbulence  and  SAS-off  cases. 
Figure  38  shows  similar  SAS  effects  for  the  rigid  body  vehicle  for  various 
turbulence  levels.  Pilot  B does  show  more  tracking  error  for  the  SAS-off 
case,  while  .pilot  C shows  a reduced  tracking  error  for  SAS  off.  Stick 
activity  data  shown  are  also  not  instructive  except  for  the  one  SAS-off 
case  for  3 fps  RMS  of  gust  velocity  which  was  described  earlier. 


Pilot  acceleration,  in  the  AR  study,  was  generated  through  the  pilot 
tracking  efforc,  the  rigid  body  motion  due  to  turbulence  and  the  acceleration 
due  to  structural  mode  excitation  by  the  turbulence  and  stick.  Figure  39 
shows  how  pilot  acceleration  is  affected  by  these  variations.  Structural 
mode  excitations  by  turbulence  are  the  dominant  contributors  for  both 
pilots.  The  pilots  really  have  no  control  over  these  effects.  The  rigid 
body  accelerations  are  considerably  different  for  the  two  pilots  for  the 
low  levels  of  turbulence,  but  are  identical  for  the  high  levels.  Although 
the  stick  activity  for  the  two  pilots  is  almost  the  same  (0.235  and 
0.225  c r ) for  3 fps  RMS  of  gust,  the  difference  in  rigid  body  accelerations 
can  be  Explained  by  examining  the  stick  spectral  density  functions.  (See 
figures  84  and  00.)  In  the  former  plot  (pilot  C)  the  power  is  concentrated 
at  about  2 rad  per  second,  while,  for  the  latter  plot  (pilot  B) , the 
dominance  is  in  the  lower  frequencies.  The  higher  frequency  inputs  by 
pilot  C result  in  attenuated  rigid  body  acceleration  effects,  while 
pilot  B's  lower  frequencies  result  in  more  acceleration  effects.  This  can 
be  seen  in  figure  40  which  presents  the  pilot  acceleration  spectral  density 
functions  for  tire  two  cases.  Pilot  B,  the  lower  plot,  shews  the  higher 
acceleration  levels  existing  at  the  lower  frequencies,  thus  giving  a 
larger  standard  de-’icu.'  n as  shown  in  figure  39.  For  the  higher  turbulence 
level  (w  = x l fps  RMSJ , the  standard  deviations  in  pilot,  acceleration 
is  about^the  same  as  seen  in  figure  40,  and  this  is  evident  in  the  similarity 
of  the  spectral  density  plots  for  tire  two  pilots,  shown  in  figure  41. 

Pilot  acceleration  data  show  little  correlation  with  pitch  stick 
activity  and  horizontal  bar  (y  R_R) , as  shown  in  figure  42.  The  arrows 
shown  on  the  curves  indicate  the  direction  of  increasing  turbulence  level. 

No  obvious  trend  exists  for  the  set  of  data,  although  some  of  the  indicated 
subsets  show  a degree  of  correlation  with  turbulence. 
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Pitch  stick  activity  trends  do  show  u|)  as  a function  of  structural  mode 
excitation.  Figure  43  shows  cr  decreases  with  gust  level  when  the  stick 
excites  the  structural  mode.  ?lie  effect  of  gust  on  the  structural  mode 
differs  between  the  pilots  with  pilot  B displaying  a reduction  with  increased 
gust  level  while  pilot  C does  the  opposite.  Both  pilots  show  a slight 
increase  in  cr  with  rigid  body  gusts.  In  the  lower  plot,  the  tracking 
error  m)  effects  are  m0re  pronounced  between  pilots  than  between  gust 
levels.  Pilot  B does  indicate  an  increase  in  cr  with  gust  increase  and 
with  gust  and  stick  excitation  of  the  structural  mode.  On  the  other  hand, 
pilot  C decreases  his  tracking  error  standard  deviation  slightly  with 
increasing  gust  level.  The  large  difference  in  tracking  error  RMS  (o'-  ) 
between  the  pilots  can  be  seen  in  the  horizontal  bar  time  histories  in 
figures  27  and  83.  In  figure  83,  it  is  evident  pilot  C allows  larger 
tracking  errors  to  exist  for  longer  times  than  does  pilot  B (figure  27). 

Analysis  of  the  pilot  bode  plots  proved  difficult  because  of  the 
variation  in  the  form  of  the  data  and  the  apparent  need  for  complex  models. 
Compared  to  the  TF  data,  the  coherence  estimates  of  the  AR  data  were  poor. 

As  described  in  appendix  I,  the  coherence  estimate  is  an  indication  of  the 
adequacy  of  a linear,  time- invariant  transfer  function  in  fitting  the 
input-output  time  histories. 

Figure  11  gives  the  coherence  function  for  a TF  run  (No.  0418)  which 
ca)i  be  compared  with  that  shown  in  figure  31  for  an  AR  run  (Nt  01A11). 

The  spectral  density  function  of  the  remnant  for  AR  run  (No.  01A11)  is 
given  in  figure  44  which  can  be  compared  to  figure  13.  The  obvious 
difference  is  the  consistently  poor  coherence  estimate  for  frequencies 
above  2.3  rad  per  second  for  the  AR  case.  This  is  typical  of  most  of  the 
AR  data,  and  examination  of  the  bode  plots  indicates  a significant  leveling 
off  or  dropping  in  amplitude  ratio  for  frequencies  above  2.5  to  3.0  rad  per 
second.  For  those  AR  cases  where  a double  lead  pilot  model  was  assumed 
and  shown  on  the  data  plots,  the  amplitude  ratio  linear  model  must  be 
restricted  to  frequencies  below  2.5  to  3.0  rad  per  second,  as  is  apparent 
from  the  plots. 

liven  in  the  lower  frequency  ranges,  a considerable  variation  in 
amplitude  ratio  form  is  indicated  by  the  bode  plots.  They  range  from 
simple  lead  (figure  66)  to  double  lead  (figure  75)  to  lead- lag  (figure  82) 
to  a notch  filter  (figure  94)  and,  finally,  to  a simple  lag  (figure  68). 

In  each  case,  the  phase  plot  does  not  match  the  amplitude  ratio  form.  The 
respective  coherence  functions  vary  as  can  be  seen  in  appendix  III.  The 
simple  lead  coherence  function  is  high  to  2 rad  per  second,  the  double 
lead  coherence  is  high  to  3.5  rad  per  second,  the  lead- lag  is  relatively 
poor  throughout,  the  notch  filter  is  poor  at  the  notch  frequency  and  above 
2 rad  per  second,  while  the  simple  lag  is  high  to  1.3  rad  per  second.  In 
general,  however,  if  the  data  points  exhibiting  low  coherence  values  are 
ignored,  the  common  tendency  is  to  a lead  only  transfer  function. 
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The  phase  plots,  in  general,  are  poor  matches  to  the  amplitude  ratio 
plots,  i-ven  where  a rather  positive  indication  of  lead  behavior  is 
evident,  as  in  figure  75,  the  phase  noticeably  exhibits  less  lead  than 
expected.  A simple  time  delay,  unfortunately,  will  not  solve  the  phase 
mismatch.  A rather  complex  time  delay  with  the  delay  time  constant,  r^}  being 
a function  of  frequency  is  required.  Run  No.  01A11,  mentioned  earlier,  is 
show  again  in  figure  45  with  a postulated  pilot  model  curve  fit  shown. 

The  frequency  response  data  can  be  matched  by  assuming  the  pilot  adopted  a 
lead  term  and  a variable  time  delay.  In  the  low- frequency  range,  he  would 
exhibit  a time  delay  of  1.0  second  which  is  removed  as  the  frequency 
increases  until,  at  3 rad  per  second,  the  time  delay  is  zero.  Figure  45 
shows  the  resulting  curve  fit.  It  can  be  expected  that  the  majority  of  the 
AR  cases  can  be  handled  by  using  a variable  time  delay  to  account  for  the 
phase  mismatch.  Without  the  benefit  of  an  automatic  curve  fit  program 
(as  was  achievable  for  the  TF  data) , such  pilot  model  construction  would 
require  lengthy  and  tedious  effort.  However,  since  the  AR  task  did  reflect 
a variation  in  pilot  rating  with  SAS  status,  another  pilot  model  construc- 
tion was  made  for  run  No.  02A2  which  was  a SAS-off  condition  with  a pilot 
rating  of  9.  Figure  46  shows  the  bode  plot  obtained  from  the  pilot  traces 
and  a pilot  model  curve  fit  based  on  a double  lead  and  a variable  time 
delay.  For  the  curve  fit,  a larger  time  delay  (1.5  seconds)  was  necessary 
than  for  the  SAS-on  condition  in  figure  45.  In  addition,  a double  lead  of 
1.0  second  was  also  necessary  for  a better  amplitude  match.  The  physical 
significar.ce  or  even  the  possibility  of  such  frequency- dependent  time 
delay  relative  to  pilot  behavior  is  not  considered.  The  pilot  remnant  PSD 
for  this  run  is  shown  in  figure  47.  No  further  attempts  were  made  to 
develop  pilot  models  for  the  various  AR  data  involving  turbulence  and 
structural  mode  variations. 

lixamination  of  the  spectral  density  functions  for  the  AR  cases  reveals 
a characteristic  which  might  be  of  use  in  correlating  pilot  ratings  with 
pilot  activity.  In  many  of  the  SAS-off  cases,  the  pitch  stick  spectral 
density  function  exhibits  a high  frequency  peak  which  is  as  large  or  larger 
than  any  existing  in  the  frequency  range  associated  with  the  tracking  task. 
Figure  36  illustrates  the  case  for  run  No.  01A11  which  is  SAS-on  with  a 
rating  of  3 and  run  No.  04A1  which  is  SAS-off  with  a rating  of  5.  Assuming 
this  is  a distinguishing  characteristic  between  good  and  bad  pilot  ratings, 
the  rest  of  the  AR  data  were  examined  for  such  traits.  For  pilot  B,  of  the 
seven  cases  with  a 5 rating,  four  checked  out  in  the  manner  previously 
described,  while  for  17  cases  with  a 3 rating,  all  checked  out.  For 
pilot  C,  all  eight  of  the  cases  rated  4 checked  out,  while  six  of  the  10 
cases  rated  2 checked  out.  Since,  in  general,  the  tracking  error  spectral 
density  function  does  not  show  a similar  behavior,  this  additional  pilot 
activity  is  similar  in  function  to  pilot  remnant  effects.  As  indicated 
earlier,  however,  the  pitch  stick  standard  deviation  does  not  correlate 
with  pilot  rating.  Further  work  in  this  area  of  spectral  density  function 
analysis  may  be  fruitful. 
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Spectral  Density  of  Pilot  Remnant  (No.  02A2) 


SECTION  VII 


DISCUSSION'  OF  PILOT  RATING  PREDICTION 


PAPER  PILOT  MODIFICATIONS 

Since  tlie  terrain- following  (TF)  simulation  did  not  corroborate  past 
investigations  relative  to  adverse  effects  of  turbulence  on  pilot  ratings, 
the  rating  algoritlim  used  in  reference  2 was  not  modified.  However,  the 
pilot  model  was  changed  to  reflect  the  characteristic  lead- lag  forn  nund 
in  the  simulation  data.  Further  modifications  were  made  to  incorporate 
the  selected  vehicle's  aerodynamic  characteristics,  its  three  primary 
longitudinal  structural  modes,  and  its  control  system  and  display  signal 
format.  The  shaping  of  the  tracking  command  signal  was  also  changed  from 
a first-order  lag  to  a second  order  so  as  to  more  closely  approximate  tire 
characteristics  of  the  simulated  7..  This  modified  paper  pilot  computer 
program  is  available  at  tire  Air  Force  Flight  Dynamics  Laboratory-,  Wright- 
Pat  terson  AFB,  Ohio. 


TERRA! N- FOLLOWING  RESUL TS 

For  the  TF  runs,  tire  pilot  lag  term  was  set  as  a constant  with  the 
pilot  gain  and  lead  time  constant  being  the  variable  parameters.  For  the 
rigid  aircraft  and  zero  turbulence  level,  tire  paper  pilot  computed  a pilot 
rating  of  1,15  with  a pilot  gain  of  4.0  inch  per  g and  a zero  lead  time 
constant.  A pilot  time  delay  of  0.42  second  was  included  ir.  tire  paper 
pilot.  These  results  can  be  compared  directly  with  TF  run  No.  1015  which, 
as  shown  in  table  IV,  has  the  same  task  conditions.  The  computed  pilot 
model  from  this  simulation  run  ..s  also  shown  in  table  IV  and  plotted  in 
figure  49.  For  convenience,  the  pilot  model  plot  is  reshown  as  figure  48 
but  with  the  addition  of  the  paper  pilot  computed  pilot  model.  The  match 
including  the  steady-state  gain  is  surprisingly  good. 

Additional  paper  pilot  runs  with  different  turbulence  levels  and  with 
structural  modes  included  showed  very  little  variation  in  pilot  rating 
and  pilot  model  parameters.  Table  III  lists  the  values  obtai’-ed  for  the 
various  combinations  evaluated.  Although  the  analytically  derived  pilot 
ratings  are  low  compared  to  those  in  the  TF  simulation,  the  same  insensi- 
tivlcy  to  turbulence  and  structural  mode  effects  is  evident. 


AIR  REFUELING  RESULTS 

The  inconsistency  in  the  pilot  model  data  obtained  in  the  air- refueling 
( ART  simulation  made  it  difficult  to  generate  paper  pilot  data  for 
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comparison.  Since  the  TF  and  AR  results  showed  little  effect  due  to 
turbulence  and  structural  modes,  only  three  runs  Qrc  shewn  in  table  III  for 
the  AR  results.  These  represent  conditions  similar  to  those  evaluated  in 
the  AR  simulation  which  resulted  in  pilot  rating  changes.  The  first  two 
AR  conditions  represented  the  change  in  SAS  command  augmentation  gains,  as 
depicted  in  figure  34.  The  primary  effect  is  in  pilot  gain  and,  not 
unexpectedly,  little  difference  in  paper  pilot  results  is  shown.  The  third 
condition  is  for  the  SAS-off  case  and  again  little  change  is  noted,  unlike 
the  simulation.  The  ratings  are  also  consistently  lower  than  those  obtained 
in  the  simulation,  again  as  was  observed  for  the  terrain  tallowing.  However, 
the  AR  paper  pilot  did  reflect  a lead  required  task  compared  to  TF  by 
reducing  the  pilot  lag  to  zero.  This  directio..  of  pilot  model  change  agrees 
with  the  simulation  data. 
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SECTION  VIII 


CONCLUSIONS 


The  examination  in  this  study  of  a large,  flexible  bomber  in  a turbulent 
environment  was  engendered  because  of  past  difficulties,  reported  by  other 
investigators,  of  correlating  increasing  pilot  rating  with  pilot  character- 
istics and  performance  parameters  when  turbulence  levels  were  increasing. 
Piloted  simulation,  conducted  in  this  study,  of  two  longitudinal  tracking 
tasks  under  high  levels  of  turbulence  revealed  little  consistent  correlation 
between  pilot  rating  and  turbulence  levels.  Pilot  acceleration,  aggravated 
by  both  turbulence  and  structural  mode  excitation,  did  not  influence  pilot 
rating.  The  pilots  participating  in  the  simulation  seemed  able  to  divorce 
their  task  rating  from  their  motion  environment.  Only  when  the  vehicle 
short  period  characteristics  were  varied  (e.g.,  SAS  off)  did  pilot  rating 
change,  an  indication  that  the  higher  frequency  effects  of  turbulence  and 
structural  mode  excitation  did  not  influence  their  primary  tracking  task. 

Pilot  performance  and  response  parameters  used  in  past  pilot  rating 
prediction  methods  did  not  correlate  with  pilot  ratings.  Pilot  acceleration, 
likewise,  did  not  provide  any  consistent  clue  to  the  piiot  ratings.  Exami- 
nation of  the  spectral  density  functions  for  the  pitch  stick  (pilot  output) 
and  horizontal  bar  (pilot  input)  revealed  a degree  of  correlation  between 
pilot  rating  and  the  pitch  stick  power  distribution  in  terms  of  frequency. 

A dominant  peak  in  the  spectral  density  function  unrelated  to  the  tracking 
task  frequency  content  seemed  to  indicate  a poorer  control  of  the  vehicle 
by  the  pilot  and,  as  a result,  a poorer  pilot  rating. 

The  pilot  models  constructed  from  the  terrain- following  (TF)  simulation 
are  considerably  different  from  those  required  for  the  air- refueling  (AR) 
simulation.  A lead- lag  model  with  a negative  lead  consistently  fitted  the 
TF  task  data  which  showed  generally  good  coherence  characteristics.  "!  e 
AR  data  varied  in  form  considerably  from  simple  leads  to  complex  lead  lags, 
to  notch  filters,  and  to  simple  lags.  In  general,  however,  if  only  data 
showing  high  coherence  are  used,  the  characteristic  linear  pilot  model  is 
of  a lead-only  transfer  function.  A frequency  dependent  time  delay  is 
required  to  provide  reasonably  good  phase  matching. 

The  results  of  the  study  indicate  that  current  flying  qualities 
specifications  should  continue  to  maintain  a separation  from  ride  qualities 
considerations  . Longitudinal  tracking  tasks  of  the  nature  evaluated  were 
rated  on  the  basis  of  the  vehicle  short- period  dynamics  and  not  on  the 
high- frequency  pilot  motion  caused  by  turbulence  and/or  structural  modes. 
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APPENDIX  I 


DATA  PROCESSING  AND  FREQUENCY  ANALYSIS  PROGRAM 


PRINCIPAL  THEORIES 

Statistical  spectral  analysis  is  a widely  recognized  tool  for  analyzing 
time-varying  data  from  a variety  of  experimental  situations.  In  this 
program,  the  emphasis  is  placed  on  cross-spectral  analysis  of  paired  time 
records  representing  an  input  time  sequence  and  an  output  time  sequence 
from  some  experimental  time- invariant  system.  Cross-spectral  analysis 
provides  a measurement  of  the  frequency  response  function  of  that  system. 

This  is  possible  even  when  die  output  is  contaminated  by  noise  disturbances 
arising,  for  example,  from  measurement  errors  or  instrumentation  variation. 

A transfer  function  can  often  be  determined  from  the  frequency  response  func- 
tion by  an  iterative  least  squares  technique.  The  principal  dieories 
(computing  formulas)  for  die  program  are  discussed  herein. 


NUMERICAL  FOURIER  TRANSFORMS 

This  program  uses  a rectangular  integration  with  end-point  corrections 
to  obtain  the  numerical  (also  called  finite  discrete)  Fourier  transform. 
Because  Fourier  transforms  are,  in  general,  complex  quantities  (i.e.,  of 
the  form  Z(jou)  = U(j«)  + jV(jw)),  die  numerical  or  finite  discrete  Fourier 
transform  may  be  represented  in  two  parts  as 


and 


N-n 

U(j<»  ) = T^  x(kT)  . cos  ko>  T 
k=n 

N-n 

V(jo>r)  = T^  x(kT)  . sin  k«T 


where 


T = sampling  period  in  seconds 
N = number  of  data  points  in  the  data  set 
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n = number  of  end  points  to  be  weighted  with  end-point  corrections 
x(kT)  - the  time  sequence  (or  correlation  sequence)  to  be  transformed 


AUTOCORRELATION  AND  CROSS- CORRELATION  FUNCTIONS 

The  following  formulas  apply  to  jointly  stationary  time  series  with 
zero  mean  values. 


Autocorrelation  of  Input 


The  input  time  series  is  denoted  bv  y(t).  The  corresponding  auto- 
correlation function  is  denoted  by  R (mT)  and  is  given  by 


m 


i 

N-m 


N-m 


Xu 


k+m 


m = 0,1,2, . . .M 


where 

m = "lag  number" 

M = maximum  number  of  lags 
T = sampling  period  in  seconds 

Note  that  this  calculation  is  for  positive  lag  numbers  only  because  R^ 
is  an  even  function,  i.e., 


R (mT)  = R (-mT) 

yy  yy 


Autocorrelation  of  Output 


'Hie  output 
autocorrelation 


time  series 
function  by 


is  denoted  by  x(t)  and  its  corresponding 
P (mT) , which  is  given  by 

XX 
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\xxW  N-m  xk  ' '^k+m 


m = 0,1,2,. . .M 


where  "M,"  and  "T"  are  as  previously  defined,  and  R is,  likewise,  an 
even  function. 


Cross-Correlation  Functions 

Botli  time  series  x(t)  and  y(tj  play  a role  in  the  calculations  of  the 
cross- correlation  functions.  The  cross-correlation  function  R (mT) 
between  output  x(t)  and  input  y(t.J  is  given  by  A’ 


R (ml)  = ~ X x.  . y, 
xy  N-m  k k+m 

K- 1 


m = 0,1,2, . . .M 


and  the  cross-correlation  function  R_  (mT)  between  input  y(t)  and 
output  x(t)  by 


R (mT) 
yx  ' 


-jl  y 

N-m  *—> 


m 0,1,2, . . .M 


These  functions  are  neither  odd  nor  even,  but  satisfy  the  relation 
(reference  8). 


R (mT)  = R (-mT) 
xy  yx  ' 


AUTOCORRiiLAT  I ON  AND  CROSS- SPbCTRAL  DLNSITY  FUNCTIONS 


Two  procedures  are  described,  each  leading  to  one-sided  estimates  of 
the  same  spectral  density  functions,  which  are  defined  for  positive  fre- 
quencies only.  The  first  procedure  is  called  the  "direct"  method.  The 
second  procedure  is  called  the  "autocorrelation"  method,  but  it  applies  to 
cross-correlation  as  well,  hither  or  both  methods  may  be  requested. 


85 


Direct  Method 

The  direct  method  starts  with  output  time  sequence  x(t)  and  input  time 
sequence  y(t)  and  obtains  the  numerical  or  finite  discrete  Fourier 
transforms  of  both  sequences;  i.e.,  it  calculates  X(jw)  and  Y(jw)  for  a 
specified  arbitrary  set  of  frequencies  o> . These  transforms  are  complex 
values  and  may  be  represented  as  real  RX(jw),  RY(jw),  and  imaginary 
IX(jw) , IY(je>)  parts  according  to 


X(j  w)  = RX(jw)  + jlX(jw) 


Y(jcoJ  = RY(jw)  + jlY(jw) 


Estimates  of  the  spectral  density  functions  are  then  obtained  as 
indicated: 


S (jw) 
xx 


[RX(jw)-jIX(ja>)]  [RX(j»)+jIX(jw)] 


S (jo)) 
xyJ 

S (jw) 
yx  J 

S (jw) 

yy 


= [RX(jw)-jIX(jw)]  [RY(jw)+jIY(iw)] 
= [RY(jw)-jIY(jw)]  [RX(jo»)+jlX(jw)] 
= [RY(jft>)-jIY(jw)]  [RY(jw)+jIY(jw)] 


By  introducing  the  notion  of  cospectral  density  C (jw)  and  quadrature- 
spectral  density  Q (jcu) , the  cross-spectral  density  terms  S (jw)  and 
S x(jft))  may  be  expressed  as  follows: 


s^yO")  = c^yO®)  + jQ^yO) 


Syx(j")  = Cyx(jW)  + jQy^  Cjw) 


it  may  be  seen  that 


C (»  = RX(j<u)*RY  (j<!))+IX(jw)  *IY  (je>) 
Cyx(jw)  = RX(j<u)*RY  (j6>)  + IX(j<u)*IY  (jo>) 

and 

Q^O)  = RX(j^)*IY  (j<y)-RY  (jw)*lX(j<y) 
A> 

QyX(j")  = RX(ja.)*IY(j£o)+RY(j£y)*IX(jct)) 

Therefore 

C (jw)  = C (jtu) 
xyw  J yx  J J 

and 

V0"1  ‘ ■V;W 


Hie  time  sequences  x(t)  and  y(t)  are  neither  even  nor  odd.  Because  the 
sequences  are  of  finite  duration,  the  end-points  introduce  integration 
errors  which  may  be  partially  compensated  by  an  end-point  correction  scheme. 
A three-point  correction  formula  known  as  Simpsons  half  rule  (reference  9) 
with  weights  5/12,  8/12,  and  -1/12  is  used  in  tire  direct  method. 


Autocorrelation  Method 


This  method  starts  with  correlation  sequences  R,  , R r>  R and  obtains 
the  numerical  or  finite  discrete  Fourier  transforms  for  ar? arbitrary 
set  of  frequencies,  on  individual  or  combined  correlation  sequences, 


as  follows: 
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V*V  - RE 


■ K [4c«  }] 

+ Im  HPq(mt)  }]  11  = M...M 

where 

PcC mt)  = 1/2  [r  (mt)  + Ryx(mt)] 

Pq( mt)  = 1/2  [Ryx(mt)  - R^Gnt)] 

As  previously  noted 

= V3"’  * V’3 

VjM)  ‘ V3"*  ' iQxy(jw) 


'Hie  autocorrelation  sequences  R (mT)  and  R (mT)  are  both  even  func 
tions  extending  to  negative  values  of* the  argumeM  mT. 

The  Fourier  transform  for  either  function,  say  R^CmT),  may  be 
expressed  as 


M 


Vn)  ■ T 2 Rx.xtn,T'cos  nT 

nl  = l 


= TR^CO)  + 2T  2 (ml) cos  mojJ 


niT=  1 


Using  a two-point  correction  for  the  end-points  with  equal  weights 
of  1/2  results  in  the  algorithm  as  given  by  equations  11  through  14  of 
reference  10.  The  same  corrections  are  used  in  this  program.  The  resulting 
expressions  are 


Rxx(0)+2  Mil  R (mT)cos  mojT 
mT=l 


+ 


R (MT)cos  MojT 


Likewise, 


R (0)+2 

y> 


M-l 

Z 

mT=l 


R^(mT)cos  mxj,T 


+ 


R (MT)cos  MujT 


Since  P (mT)  is  an  even  function,  its  Fourier  transform  may  be  obtained 

in  the  same  manner  as  for  R and  R previously  noted,  i.c., 

xx  yyJ  / > > 


P (0)+2  ^ P (mT)cos  mo^T 
mT=l 


+ 


as 


Since  P^(mT)  is  an  odd  function,  its  Fourier  transform  may  be  expressed 


= 2 


M-l 

0+2  £ P (mT)sin  m«T 

. mT=l  q % 


+ P^(MT)sin  Mc^T 


Variability  of  Spectral  Density  Estimates 

This  discussion  applies  to  the  autocorrelation  method  only. 

Variability  of  the  spectral  density  calculations  is  said  to  be  chi- 
squared  distributed  with  degrees  of  freedom  2N/M  when  S and  S are 

suitably  normalized.  This  drops  to  N/M  at  frequencies  0 an^oo  = vr/T. 
The  parameters  N,  M,  and  T are,  respectively,  the  number  of  time  data 
points,  the  number  of  lags  used  to  calculate  the  autocorrelation  sequence, 
and  the  sampling  period  in  seconds.  The  required  conditions  for  a chi- 
squared  variate  are  that  x(t)  and  y(t)  be  Gaussian  and  their  ideal  spectral 
density  functions  and  S be  reasonably  constant  in  the  interval 


l 


U<0), 


These  considerations  rest  with  the  user  to  interpret  for  his  own  use. 
Refer  to  Part  1,  section  14,  of  reference  11  for  more  detail. 

This  program  was  designed  so  that  the  ratio  of  N/M  would  not  be  less 
than  2 whenever  the  autocorrelation  method  is  specified.  This  is  the 
recommended  minimum  for  that  ratio  and  is  equivalent  to  two  degrees  of 
freedom  in  the  chi-square  sense.  The  program  calculates  the  value  of  M 
from  the  punched  card  entry  for  DF,  the  frequency  increment,  according  to 


M = 


I 


7T 

T*DF 


M < 500 


I' 


90 


Since  the  program  requires 


2 


then,  substituting  from  M - and  rearranging  gives  the  required  range  for 
DF,  i.e., 


DF 


27 T 
- N*T 


N < 10,000 


If  DF<27r/N*T  the  program  will  increase  its  punched  card  value  25  percent, 
i.e. , 


Then,  if 


DF’ 


2 7T 

N*T 


tiie  program  will  substitute  the  direct  method  for  the  autocorrelation  method 
and  continue. 

Very  often,  a choice  in  sampling  rate  fg  for  the  time  sequences  x(t) 
and  y(t)  may  be  exercised  before  the  data  are  digitized.  For  data  sampled 
several  factors  above  the  bandwidth  of  the  system  producing  the  data,  the 
program  offers  a simple  reduction  feature.  However,  one  must  be  careful 
not  to  alias  the  data  by  sample  reduction.  Aliasing  is  said  to  occur  when- 
ever the  sampling  frequency  is  less  than  twice  the  highest  frequency  of 
interest.  For  example,  consider  a low  bandwidth  s, stem  bandlimitea  at 
10  radians  per  second.  The  slowest  sampling  rate  that  avoids  aliasing  is 


f = 
s 


10 

IT 


3 samples  per  second 


SMOOTHING  TO  IMPROVE  SPECTRAL  ESTIMATES 


in  estimating  spectral  densities,  an  ideal  representation  for  the  time 
function  x(tj  is  a record  of  infinite  length.  In  practical  situations, 
only  short  intervals  of  time  are  available  for  analysis.  Also,  e may 
require  shorter  intervals  of  time  for  other  considerations  than  are 
actually  available.  A mathematical  means  which  literally  accounts  for 
"turning  the  time  function  on  and  off"  is  to  use  a spectral  filter  or 
window.  A so-called  Hanning  window  is  used  in  this  program  (reference  11). 
This  filter  smooths  the  spectral  density  estimates  across  three  spectral 
points. 


COHERENCE  ESTIMATES 

Coherence  is  a measure  of  the  degree  to  which  two  time  series  x(t)  and 
y(t)  are  related  by  a linear  model  of  the  type  shown  here. 


(Time- Invariant)  N(t) 


where  N(tj  is  extraneous  noise.  The  two  time  series  x(t)  and  y(t)  are 
supposed  to  be  jointly  "weakly"  station'3'*''  (i.e.,  the  values  of  their  first 
and  second  moments  are  independent  of  ti^  choice  of  the  time  origin) . The 
noise  disturbance,  N(t),  must  be  uncorrelated  with  the  input.  For  classical 
correlation  analysis,  the  linear  model,  G^(jw)  may  be  viewed  on  a con- 
ceptual basis  only.  In  any  event,  cross -spectral  analysis  seeks  to 
determine  the  degree  to  which  a linear  relationship  prevails  between  two 
time  series  and  to  characterize  the  linear  model  (or  filter)  which  may  be 
thought  of  as  providing  this  relationship.  Then,  in  ter^s  of  the  spectral 
density  functions,  the  coherence  functions,  denoted  by  R (a/),  is  given  by 


R2(u>)  = 


C2  (jw)  + Q2  (j  (J) 
xy  J x xy  J 

S (jta)S  (j<J) 
xx  J yy  J 


As  long  as  the  linear  model  holds,  the  coherence  function  can  also  be 
expressed  as 


R2(M  = 


1 


S 

nn 


S 

xx 


CM 

(M 


where  S (jM  is  the  spectral  density  function  of  the  output  disturbance, 
N(t).  Tne  total'  spectral  density  of  the  output  can  be  written  in  terms  of 
snn(j^  and  Sx*x*(;ja>)  as 


S 

xx 


CM 


X*X 


(jw)  + SnnM 


where 


S^CM  ~ the  spectral  density  of  x*(t) 


Thus,  it  is  apparent  that  coherence  always  assumes  values  between 
zero  and  one.  In  the  absence  of  an  output  disturbance  (i.e.,  if 
^CM  = 0 at  some  frequency  uj  then  R (u  ) = 1,  and  x(t)  and  y(t)  are 
"perfectly  coherent"  (at  that  triquency.  The^opposite  condition  is  where 
(Wj<.)  = 0,  then  x(t)  and  y(t)  are  "perfectly  incoherent.") 

It  should  be  noted  there  is  no  theoretical  significance  to  the 
preceding  calculations  of  coherence  when  the  direct  method  is  specified 
since  the  equation  results  in  values  of  unity  for  all  frequencies. 


GAIN  AND  PHASE  ESTIMATEES 


Gain  and  phase  estimates  are  obtained  from  the  linear  model  as  the 
ratio  of  X(jw)/Y(joj),  i.e., 


gaCM  = XCM/YCM 
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where  G.  (jo>)  may  be  regarded  as  the  frequency  response  function  relating  the 
input  ytt)  to  tire  output  x(t).  It  may  also  be  regarded  as  the  transfer 
function.  If  we  multiply  both  numerator  and  denominator  by  the  conjugate 
of  Y(jo»),  denoted  by  Y(jw),  then 


GA(JWJ  Y(jcJ)*  Y(jo>) 


Y(jw)  * X(jw)  = Crx(j^  + jQ^CW) 


Y(jw)  * Y(jcJ)  = S^Cjw) 


c CjcJ)  + jO  ( 
V*  ■ y S (j 

yy 


Therefore,  the  gain  of  this  transfer  function  is  computed  as 


GA(jw) 


S (j(j) 

yy  " 


or  since 


C = C and  0 = -Q 

yx  xy  yx  Nxy 


APPENDIX  II 


DATA  ANALYSIS  RESULTS  FOR  TERRAIN- FOLLOWING  SIMULATION  RUNS 


The  extensiveness  of  the  data  presented  in  this  appendix  is  to  provide 
a complete  data  base  for  the  simulation  tests. 

A listing  of  the  terrain- following  (TF)  runs  and  the  significant  data 
analysis  results  are  given  in  table  IV.  Figures  49  through  58  present  the 
data  output  plots  obtained  from  the  computer  data  processing  of  die  listed 
runs.  Each  figure  shows  die  pitch  stick  and  horizontal  bar  (ADI  flight 
director  horizontal  needle)  time  histories;  the  autocorrelation  functions; 
the  spectral  density  functions;  die  cospectra,  quadrature  spectra,  and 
coherence  functions;  and  finally  die  bode  plots  wid'  die  computed  curve  fit 
transfer  functions.  In  many  instances,  it  can  be  seen  diat  a better 
amplitude  ratio  fit  can  be  obtained  if  die  curve  fit  is  lowered  in  the 
amplitude  plots  without  any  phasing  changes.  This  corresponds  to  just 
lowering  the  dc  level  of  the  transfer  function.  The  values  listed  in 
table  IV  reflect  these  changes. 

Figure  59  presents  data  generated  using  the  pilot  model  as  a means  of 
closing  die  g control  loop  for  the  TF  task.  The  specific  model  used  is 
that  listed  under  run  No.  0418  in  table  IV.  Computer  runs  were  with  diis 
pilot  model  and  also  with  an  automatic  TF  system.  The  graphs  indicate  die 
pilot  model  performed  quite  satisfactorily. 
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DATA  ANALYSIS 
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1015  MANUAL  terrain  following  analysis. 
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Figure  49.  TF  Simulation  Run  No.  1015  Data  (Cont) 
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Figure  49.  TF  Simulation  Run  No.  1015  Data  (Cont) 
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Figure  50.  TF  Simulation  Run  No.  0416  Data  (Cont) 
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0417  MANUAL  TERRAIN  FOLLOWING  ANALYSIS 
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Figure  51.  TF  Simulation  Run  No.  0417  Data 
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Figure  52.  TF  Simulation  Run  No.  0418  Data  (Cont) 
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Figure  54.  TF  Simulation  Run  No.  0248  Data  (Cont) 
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Figure  54.  TF  Simulation  Run  No.  0248  Data  (Cont) 
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Figure  54.  TF  Simulation  Run  No.  0248  Data  (Cont) 
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Figure  54.  TF  Simulation  Run  No.  0248  Data  (Cont) 
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Figure  56.  TF  Simulation  Run  No.  0254  Data  (Gont) 


HORIZONTAL  BAR  AND  PITCH  STICK  - IN.  X G PITCH  STICK  AND  HORIZONTAL  BAR 


025H  MWJUAL  TERRAIN  fOLLOHING  ANALYSIS. 
CROSS  CORRELAT ION  fUHCTIONS 


° .02 


LAG  NUICER 


Figure  56.  TF  Simulation  Run  No.  0254  Data  (Cont) 


145 


/SEC 


COHERENCE  QUADRATURE  SPECTRA  COSPECTRA 


0204  MANUAL  TERRAIN  fOLLOWlNG  W YSIS. 


AUTO-CORRELAT ION  KE1H0D 


.12 


.08 


.04 


1 

| 

1 

1 

II 

\mm 

i 

■ 

i 

1 

8 

1 

I 

II 

II  I 

i 

! 

■ 

1 

■ 

1 

1 

1 

II 

■ 

■ 

I 

i 

■ 

i 

■ 

■ 

1 

1 

II 

i— B 

u 

■ 

■ 

■ 

i 

i 

i 

i 

5 

B 

i 

i 

■ 

| 

■ 

■ 

H 

| 

I 

II 

IB 

u 

B 

■ 

■ 

i 

i 

i 

i 

PI jW 

■ 

■ 

■ 

i 

i 

I 

i 

1! 

■ 

■ 

1 

| 

I 

II 

■S 

m 

■ 

■ 

E 

i 

i 

1 

i 

ma 

■ 

E 

■ 

i 

i 

g 

1 

B 

■ 

8 

1 

1 

8 

m 

B 

■ 

■ 

i 

i 

i 

i 

■ 

■ 

i 

i 

8 

1 

■ 

1 

II 

IS 

u 

■ 

■ 

i 

i 

i 

i 

■ 

■ 

■ 

i 

i 

i 

1.X10 


-01 


,♦00 


1 .X10  ””  1.X10 

FREQUENCY  (RAO/SECl 


♦ 01 


.100 


1.0 


Figure  56.  TF  Simulation  Run  No.  0254  Data  (Cont) 


145 


I . X 1 0 


♦ 02 


♦ 02 


'02 


..  


I 

/*> 


1 

J 


PHASE  - DEG  GAIN 


0254  MANUAL  TERRAIN  FOLLOWING  ANALYSIS. 


to 

O 


15 


Figure  56.  TF  Simulation  Run  No.  0254  Data  (Concl) 


146 


fa*!  1 T mi' 





HORIZONTAL  BAR  - G2  P,TCH  STICK 


148 


1 


0256  MANUAL  TERRAIN  FOLLOWING  ANALYSIS. 
SPECTRAL  0ENS1IY  FUNCTIONS 


Figure  57.  TF  Simulation  Run  No.  0256  Data  (Cont) 
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Figure  o9 . Comparative  TF  Runs  With  Pilot  Model  and  Automatic  TF  System  (Concluded) 
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DATA  ANALYSIS  RliSULTS  FOR  AIR  REFUELING  SIMULATION  RUNS 


As  for  the  terrain- following  (TF)  data,  the  extensiveness  of  the  data 
in  this  appendix  is  to  provide  a complete  data  base  for  the  simulation 
tests. 
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Identification  and  data  results  for  the  air  refueling  (AR)  simulation 
runs  are  listed  in  table  V.  Figures  60  through  99  show  the  data  plots 
obtained  from  the  computer  processing  of  the  simulation  data.  Each  figure 
shows  the  pitch  stick  and  horizontal  bar  time  histories;  the  autocorrelation 
and  cross- correlation  functions;  the  spectral  density  functions;  the 
cospectra,  quadrature  spectra,  and  coherence  functions;  and,  finally,  the 
bode  plots  for  the  pilot  transfer  functions.  As  indicated  earlier,  no 
automatic  curve  fitting  was  achieved  and,  therefore,  none  are  shown  on  the 
bode  plots. 
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Figure  60.  AR  Sinrulation  Run  No.  01A12  Data  (Cont) 
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Figure  61.  AR  Simulation  Run  No.  01A13  Data  (Cont) 
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Figure  61.  AR  Simulation  Run  No.  01A13  Data  (Cont) 
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Figure  62.  AR  Simulation  Run  No.  01A14  Data  (Cont) 
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Figure  63.  AR  Simulation  Run  No.  01A15  Data  (Cont) 
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Figure  63.  AR  Simulation  Run  No.  01A15  Data  (Cont) 
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Figure  64.  AR  Simulation  Run  No.  01A16  Data 


190 


PITCH  STICK  AND  HORIZONTAL  BAR 

HORIZONTAL  BAR  AND  PITCH  STICK  - IN.  X DEG  IN.  X DEG 


^MBBBBBBBBBBBI 

■KBBBBBBBBBBBBBflBBBflBflBBBBBBBBBBBBBB 

■■■■■■■—■■■■■■■■■■■■■■■■■■■■■■■■■■■I 

■■  i r mmmmmmmm 

■■■phrspbI 

mmmmmmmmm 

■BBBBBBBBH 

. a 

IBB 

mm 

ISSbbmbSb 

pbSSbbbbbSSSS| 

Bbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbi 

SMBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB^HBl 
^SmmBgflBBBBBBflBBBBg^^^^^^^^^^^gj 

IflBBUHBBBBBBBBBBflBflBBBBBBflHBflBflBm 
■BBBBBflflBBBBBBBBBflBBBBBBBBBBBBPBBBBB 
■BBBBBSgggggBBBBBBBBBBBBBBBBIgggBBBBMl 

■■■■■BHaBBBHiHaaHi ■ 

Bfl 

BBfll 

^^^^■iBBBBBBBBBSiBBBBBBflBBBBBBBflBBB| 
^^^^■BBBBBBBBBBBBBBBBBBBBBBBBBBBBBl 


'IP 


COHERENCE  QUADRATURE  SPECTRA  COSPECTRA 


194 


GAIN 


HORIZONTAL  BAR  - DEG  PITCH  STICK 


HORIZONTAL  BAR  - DEG 


01A17  REFUEL  ANALYSIS 
AUTO  CORRELATION  FUNCTIONS 


■■■■■■■■■■■■■■I 

^^^^■BBBflBBBBBBBBBBBBBBBBBBflBBBBBBBBBBBB 

^^^Mbbbbbbbbbbbbbbbbbbbbbbbb ■■■■■■■■■«■! 

^^^HflBBflflflBflBBflflBBflBflBBflBflBBBBMHHHMHB 

^^^^■llflflBBBBBBBBBBBBBBBBBBBflBBBflBBBflflBBB 

^^^■BflBBBBBBBBBBflBBBBBBflBBBBBflflflBBBBBBB 
klBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB 
■KBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBB 
BBgBggBBQBBBBBBBBBBBBBBBBBBBBBB 
■■■■■iilBIBI  BBBBSSIBBBBieSSIBflBi 

BBBflSBflBfilflBIfiBBlifliBBBBB 

bbSiiiSbbbSSbbbSSSSiiS 

bssbssBh 

HH 

!9!nnn^BB!ln! 

SSSiSiSSSSiSiiiSSSiSiSSSiSiiiiiii 

flIIIIIIIIHIIIIIIIIII 
■KIHIIHIIIHHIIllI 
HHHHIbbi  isHkiummi 

uni  ii  lauiiiiniini 


LAG  NUMBER 


LAG  NUMBER 


Figure  65.  AR  Simulation  Run  No.  01A17  Data  (Cont) 
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Figure  67.  AR  Simulation  Run  No.  01A23  Data  (Cont) 
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Figure  68.  AR  Simulation  Run  No.  01A24  Data  (Cont) 
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Figure  68.  AR  Simulation  Run  No.  01A24  Data  (Cont) 
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Figure  69.  AR  Simulation  Run  No.  01A33  Data 
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Figure  70.  AR  Simulation  Run  No.  01A34  Data  (Cont) 


228 


01A3M  REHXL  ANALYSIS 


10 


Figure  70.  AR  Simulation  Run  No.  01A34  Data  (Concl) 
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Figure  71.  AR  Smilation  Pun  No.  01A35  Data  (Cont) 
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Figure  72.  AR  Simulation  Fun  No.  01A38  Data  (Cont) 
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Figure  72.  AR  Simulation  Run  No.  01A38  Data  (Cont) 
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Figure  73.  AR  Simulation  Run  No.  01A39  Data 
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Figure  74.  AR  Simulation  Run  No.  01A40  Data  (Cont) 
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Figure  74.  AR  Simulation  Run  No.  01A40  (Concl) 

255 


—a.. ....—,  i i 


TIME  (SECONDS) 


Figure  75 


t 


. AR  Simulation  Run  No,  01A3  Data 


256 


■ .A. 


01  A3  REFUEL  ANALYSIS 
AUTO  CORRELATION  FUNCTIONS 
.10 


it 

5* 

r 

F 

?■ 

( 

r 

*1 

!»■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■! 
MBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBl 
IIBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBBH 

llllBBIBBBBBBBBBBflBBIBBBBBBBBBIHiiB 
■iSilBaBBaBHlBBBBlWBlBIBlMBBlBB! 

■■BBBBMIWfHBBBBBBBBBBBBBBBBBBI 


LAG  NUMBER 


HBmbbmbmmhMH 

liiiiiiiiiimiiiiiiiiiiiniiiiniii 

ikriiiiiiiiirrrrrhhrrrhiirrhrrrihir 

■ikiiiiiiiiiiiikisbiibirsir&ihiiiBi 

RRRnRRRRRBaRRRSRRRRKRRRRRRRRRRHRRRg 

RRRKRBRRRRRRRRRRRRRRRRRRRRRRRRRRRRR 
RRRRBRRRRRRRRRRRRRRRRRRRRRRRRRRRRRR 
— RRRRRRRIRRIRRlBlwim 


LAG  NUMBER 


Figure  75.  AR  Simulation  Run  No.  01A3  Data  (Cont) 
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Figure  75.  AR  Simulation  Run  No.  01A3  Data  (Cont) 
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Figure  76.  AR  Simulation  Run  No.  01A10  Data  (Cont) 
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Figure  77.  AR  Simulation  Run  No.  04A1  Data  (Cont) 
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Figure  78.  AR  Simulation  Run  No.  0?A2  Data  (Cont) 
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Figure  79.  AR  Simulation  Run  No,  01A19  Data  (Cont) 
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Figure  80.  AR  Simulation  Run  No.  01A20  Data  (Cont) 
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Figure  81.  AR  Simulation  Run  No.  01A31  Data  (Cont) 
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297 


HORIZONTAL  BAR  - DEG  PITCH  STICK 


' * / 


m 

HI 


IllHHHI 

IllHHHI 

■■■■■■■I 


ini 


Hill |g|  m ■ ■ ■ Va UU  m m ■■  fUH #s •»  arm- mi 


iiiii 

■■■■■ 

«■■■■ 

huh 

■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■ 

■■■■■■■■■■■■■■■■■■■■■»■■■■■■■ 

■■■■■■■■■■■■■■■■■■■■■■■■■■■■■■; 

HUHHHHH 

HHUHHHH 

HHMHH 

BSBwSSSS 

■■■■■ 

■■■■■ 

niai 

■■■■■ 

■■■■■ 

■■■■■ 

■■■■■ 

■■■■■ 

■■■!■ 

■■■■■ 

■■■■■ 

■■■■■ 

■■■■■ 

■■■■■ 

■■■■■ 

■■■■■ 

■■■■■ 

iiiii 

■■■■■ 

■■■■■ 

■■■■■ 

■■■■■ 

■■■■■ 

mmmmm 

88888 

■■■■5 

BESS* 

HU 

iaa 

■■■■■■■■■■■■■■■■■■■■■■■■■■■ 

8I8888888888IB8SS88SE888888I 

mi 


iBssssa 
!■■■■■■ 
mmmmSm 


SH 


01 A3£  fiCFUEL 


■■III 


01 A32  REFUEL  ANALYSIS 


t 


M 

Si 

;1 


1 


? 


i 


t 

Figure  82.  AR  Simulation  Run  No.  01A32  D ,ta  (Concl) 


303 


Jf  <A-V 


j 

1 


f 

i 


i 


I 


08A49  REFUEL  ANALYSIS. 
SPECTRAL  DENSITY  FUNCTIONS 


AUTO-CORRELATION  METHOD 


>1 


i 


4 


■s 


i 


1 

i 
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Figure  84.  AR  Simulation  Run  No,  01A50  Data  (Cont) 
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Figure  84.  AR  Simulation  Run  No.  01A50  Data  (Cont) 
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Figure  85.  AR  Simulation  Run  No.  01A51  Data  (Cont) 
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Figure  86.  AR  Simulation  Run  No.  01A52  Data  (Cont) 
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Figure  96.  AR  Simulation  Run  No.  01A54  Data  (Cont) 
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Figure  96.  AR  Simulation  Run  No.  01A54  Data  (Cont) 
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Figure  96.  AR  Simulation  Run  No.  01A54  Data  (Concl) 
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Figure  99,  AR  Simulation  Run  No.  01A62  Data  (Cont) 
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Figure  99.  AR  Simulation  Run  No.  01A62  Data  (Cont) 
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Figure  99,  AR  Simulation  Run  No.  01A62  Data  (Concl) 
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